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SECTION 1.0 
INTRODUCTION 


A filter is desired to separate the geoid undulation 
information from measurement noise in the GEOS- 3 altimeter 
data. Obviously, for any separation to be possible at all, 
the geoid undulation and the noise must possess some distinct 
characteristics. The synthesis of the filter and the result- 
ing filter performance depend on the available information 
about these characteristics. Frequently the availability 
and reliability of this information are the weakest link in 
the data filtering process. 

In this analysis, it is assumed the geoid undulation 
and the altimeter noise are uncorrelated random processes 
with known second order statistical properties (i.e., means 
and correlations) . The modeling of the signal (geoid undula- 
tion) as a random process is comparatively recent and its 
conceptual basis was clarified in the fundamental work of 
Shannon (Reference 1) . The restriction to second order 
statistics arises out of necessity because of a lack of 
knowledge of the underlying physical processes. It is 
expected the analysis of GEOS-3 altimeter data will yield 
improved information about the spatial auto- correlation of 
the geoid undulation. Before such analysis has been com- 
pleted, however, the TASG model (Reference 2) seems to be 
the best available, and will be used. The altimeter noise 
shall be assumed correlated exponentially in time. At low 
altimeter sampling rates when the sampling period is long 
in comparison with the noise correlation time, the noise 
may be considered uncorrelated from sample to sample. 
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1 
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As long as higher order statistical information is not 
available, an optimum linear filter may be considered the 
optimum filter. The proposed filter will consist of a set 
of computational algorithms easily implemented on a digital 
computer. The derivation of the optimum filter algorithms 
is given in subsequent sections. A realistic numerical 
example showing the convergence characteristics of the 
filter and accuracies to be expected from the filtering 
operation is given in Section 5.2. 
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SECTION 2.0 

THE OPTIMUM LINEAR FILTER 

The GEOS-3 altimeter measurements, after the orbit 
information and other known corrections have been subtracted 
may be written as 


7i = hi + n^ (i=0,l,2,3. . .N-1) 


C2.1) 


where 

h = geoid undulation 

n = measurement noise 

i = index for measurement time. 

The objective of the filter is to determine an optimum esti- 
mate of the geoid undulation h*(t) from the N measurements 
y^ (i = 0,, 1 , 2 , 3 , . . ,N- 1) . A general linear estimate may be 
written as 

N-1 

hCt) = ^ upt) y. (2.2) 

i=0 

with an estimate error 

N-1 

e(t) = h(t) - h{t) = ^ to^(t) - hCt) (2,3) 

i=0 


3 


The weighting factors w^(t) are the filter coefficients which 
define the filter. The optimum filter is obtained if one 
may impose the following condition on the weighting factors 

N-1 

E(eCt)y^) =^oj.(t3 E(y.y^) -E(h(t)y^)=0 (2.4) 

i=0 

(il=0,l,2. . .N-1) 

Equation (2.4) states there exists no correlation between the 
estimate error and the measurements, and therefore there is 
no way (within second order statistics) of extracting more 
information out of the measurements. Since the measurement 
autocorrelation matrix E(y-yj) may generally be taken to be 
positive definite. Equation (2.4) may be inverted to obtain 
the optimum weighting factors as, 

T "1 

[u*Ct)l = E([y][y]‘) E(h(t)[y]) (2.53 

Nxl NxN Nxl 

where 

[ui C^)j ” ^^0 ^2 ••• ^N- 1 

[yf ' ^2 ••• ^N-d 

The meaning of the weighting factor as given by Equation (2.5) 
is obvious and is intuitively reasonable. Notice that these 
optimum weighting factors depend only on the correlation func- 
tions and do not presume the existence of a functional relation 
between the signal and the measurement such as Equation (2.1). 
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Under the assumption that the geoid undulation and the 
altimeter noise are uncorrelated, one has the further 
s implif ication 


Eay][y]’^) = E([h][h)h + E([n](n]'^) (2.6) 

E(h(t)[y]) = E(h(t)[h]) (2.7) 

[h] - [Eq ii2 ... 

Equations (2.2), (2.5), (2.6) and (2.7) show that formally 

the optimum filter can be constructed when the auto- 
correlation function of the geoid undulation and the altimeter 
noise are given. However, the algorithm given by these equa- 
tions as they stand involves a great deal of computation and 

is impractical for processing a large amount of data. To 
overcome this difficulty one needs to postulate simple 
structures for the underlying signal and noise processes. 

Since all one knows about these processes are second order 
statistics, it is only important that the postulated pro- 
cesses have the same second order statistics. Wiener 
(Reference 3) and Kalman (Reference 4) showed us that if 
such processes may be generated as the outputs of finite- 
dimensional linear dynamic systems (e.g,, ordinary dif- 
ferential or difference equations or the corresponding 
transfer functions) driven by white noises, computationally 
efficient algorithms may be readily devised to implement 
the optimum filter. In other words, one looks for finite- 
order Markov process representation of the signal and the 
noise. Being Markovian means only a small number of parame- 
ters need be kept, resulting in vastly reduced computational 
and memory requirements . 


If the correlation functions are stationary, the Levinson re 
cursion may be used to lessen the computational requirement 
(Ref. 3) . 


From the above discussion one sees that given the geoid 
undulation correlation function, the key to the construction 
of the optimum filter for the GEOS- 3 altimeter data is the 
modeling of the geoid undulation as the output of a linear 
ordinary differential equation system driven by white noise. 




SECTION 3.0 

MARKOV MODELS FOR GEO ID UNDULATION 
AND ALTIMETER NOISE 

3.1 WHITE-NOISE-DRIVEN DIFFERENTIAL EQUATION MODEL FOR 
TASC GEOID UNDULATION CORRELATION FUNCTION 

The Analytical Sciences Corporation (TASC) introduced 
the following model for the correlation of the geoid undula- 
tion with surface distance (Ref. 2). 


E(h(5-+d) h(d)) = a\i + XZ + 




(il>0) 


(3.1) 


Since GEOS- 3 is in a circular orbit with "constant” ground 
speed, one may convert this to a correlation in time as 

E(h(x+t) h(t)) = a2(i+3T+ — T^) (3.2) 

3 

(t>0) 

First of all, notice this correlation function is stationary; 
i.e., it depends on the time separation t only. This means 
the correlation function may be the output of a differential 
equation with constant coefficients. Secondly, the correlation 
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function depends on three linearly independent functions of 
time, implying it needs a 3rd order system to describe it. 
Thirdly, the forms of these linearly independent functions 
e ^ , Te ^ and x e ^ are reminiscent of a linear differen- 
tial equation with triple roots -3, or the geoid undulation 
may possibly be modeled by the system 

+ 8)^ h = w (3.3J 

where w is a zero-mean white noise with correlation function 

E(w(t)w(x)) = q''6(t-x) (3.4) 

It will be shown that indeed the above statement is correct. 

With knowledge of the differential Equation (3.3) or the 
corresponding transfer function 

H(s) = l/(s+8)^ 

the classical procedure would be to design a Wiener Filter 
implemented by hardware. The Wiener Filter tacitly assumes 
the filtering process has reached a "statistical steady state" 
and may involve some errors for short GEOS- 3 passes. Since 
the present filter is to be implemented on a digital computer, 
no such sacrifices are necessary and we shall adopt the 
Kalman formulation in state space. To proceed, let us 
define a convenient set* of state variables as 

'S 

Different choices may be made. 
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-Xi" 


+ 30^ h 

^2 

= 

(^ - 30 h 

-^3-J 


- h 


which satisfy the state equation 



"3 0 o' 


X 

1 


1 

1-3 0 


^2 

+ 

0 

0 1-3 


_^3_ 


0 


(3.60 


It may be shown by substitution that these state variables 
have the state transition matrix 


[Kt,O0] 




(3.70 


The general solution of Equation (3.60 is 


[x(t01 = [nt,O0] [x(O0] + / e 


i 


^-3 (t-rO 


t-T 

.jCt-T) 


w(xOdT (3.80 


From Equation (3.80 and the fact that future input is uncor- 
related with past state, one obtains the correlation matrix 

E([x(tO] [x(O0]'^)= [Kt,O0] E([x(O0] [xCOOl^O 

q 


In particular. 


ECXjCtJXjCO)) = {ECx3(0)2)+t E(X2(0)X3(0)) 


+ ECXjCOXjCO))} 


(3.9) 


Since the above expression is nothing but the geoid undulation 
correlation function, by comparing it with Equation (3.2), one 
must have 

E(x3(0)^) = 0^ 

23^ 2 

E(Xt (0)x,(0))= 0 (3.10) 

-^ 3 

E(x2(0)x3(0)) = 30^ 

Alternatively, one may obtain the covariance matrix of the 
state variables from Equation (3.8) as 

E([xCt)l [x(t)jh= t'Kt.O)] Eax(0)][x(0)f)[(Kt,0)]'^+[Q(t,0)] 

(3.11) 

where 
[Q(t,0)]4 

(3.12) 

is the increment in the covariance matrix due to white-noise 
excitation. The evaluation of the integral gives us the 
elements of [Q] as 


^11 X. 


^12 ^2 2 


q'/o e-^e(t-t) 


^13 ^23 ^^33 




t-T 

.|(t-x)2 


[i! t-xi ] 


dx 
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N 


I 


r J 


1 * 


1 


^11 

= Y(l-e 

-2Bt^ 


^12 

= aqil" 

yte-2^^ 


^13 

= °‘^12" 

0.5 Yt^e 

-2gt 

^22 

= 2qi3 



^23 

= 3aqj^3 

-0.5 Yt^ 

^-2Bt 

^33 

2aq23 

-0.25 yt 

4^-2St 

(7 

= qV23 

, a = 

1/2S) 


(3.13) 


which consist of a constant part and a damped transient part. 

For the system to be stationary, the white noise excited time- 
dependent transient part must be balanced by the free transients 
of the system. It is straightforward to show from Equation 
(3.11) that this is indeed so, provided q = 3 — and then 

one has , 

E([x(t)Hx(t)]''') [QCt.O)] (3.14) 


Notice this result is in agreement with those in Equation (3.10). 

It has been shown that the state variable of the 
linear system Equations (3.5) and (3.6), with the a priori 
covariance, Equation (3.14), models the geoid undulation with 
the correlation function. Equation (3.2). One may give the 
following loose but physically meaningful interpretation to 
the state variables defined in Equation (3.5), and to 
Equation (3.6) governing these variables. x^, ^2 ^1 

related to the geoid undulation, its ’’slope'* and "curvature", 
respectively. Equation (3.6) says the rate of change of the 
curvature of geoid undulation, although having similar mean 
square values, is uncorrelated from one location to another. 
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3.2 


ALTIMETER NOISE MODEL 


Two models of the altimeter noise sequence n^^ are 
considered. For sampling rates at 10 samples/sec or lower, 
the altimeter noise will be considered uncorrelated from 
sample to sample; i.e.. 


ECn.n.J 


= R. 6. . 
1 ij 



j=i 


(3.15) 


For higher sampling rates the noise sequence will be considered 
to be correlated exponentially with some correlation time 1/b; 
i.e., 

-b(t.-t .) 

E(n. n^) = R. e ^ J (i>j) (3.16) 


The exponentially correlated noise may be modeled, similar to 
Xj in the preceding section, by the following Markov sequence. 


n. = e 


-b (t . - t . , ) 
^ 1 1 - 1 '^ 


n. , + u. 
1-1 1-1 


E(u^uj^) = R^(l-e 




) 


£ > k 


(3.17) 


When the correlation time approaches zero; i.e., b-><», the 
exponentially correlated noise approaches the white noise. 
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SECTION 4.0 

OPTIMUM ZERO- LAG RECURSIVE FILTER (KALMAN) AT LOW 
ALTIMETER SAMPLING RATES (10 s.^mples/sec) 

The measurements, the state variables governing the 
geoid undulation and the altimeter noise are described by- 
equations in the preceding sections in the standard forms 
for an immediate construction of the celebrated Kalman filter 
(Reference 4). When no time lag exists; i.e., estimates of 
geoid undulation at any time t are to be based on past meas- 
urements only, this filter is a recursive realization of 
the optimum filter. Equation (2.5). When all data are 
considered, the zero- lag filter is not optimum because it 
does not make use of later measurements. However, it is 
known that an optimum "smoother”, making use of all the data, 
may be constructed by properly combining a forward and back- 
ward pass of the data through a zero-lag filter (Reference 5) , 
and this will be discussed in Section 5. 


A simple derivation of the zero-lag recursive filter 
for altimeter data is given in the following for the benefit 
of those who may not be familiar with Kalman filter and for 
the physical insight the derivation provides. In the deriva- 
tion we shall make constant use of the following: 


1. White noise is "purely random”, or unpredictable. 

2. The optimum estimator. Equation (2.5), is reduced 
to the following for a single measurement. 



optimum estimate of signal 


Cross-correlation of signal 
and measurement 

Auto-correlation of measurement 


* (measurement) 


(4.0) 
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4.1 A PRIORI ESTIMATE, INTERPOLATION BETWEEN AND EXTRAPO- 
LATION BEYOND MEASUREMENTS 

Prior to any measurements, the geoid undulation and 
the other related state variables are considered zero-mean 
random processes with covariance matrix Equation (3.14). 
Therefore one has 


[x(0")] = [0] 


(4.1) 


with 


E([x(0')-x(0)][xC0:)-x(0)]') = [P(0")] =Y 


la a 

2 3 

a 2a^ 3a'^ 

3cc^ 6a^ 


C4.2) 


The geoid undulation state at two time instants, t and 


t^£t are related by 


- 1 


<^1 


1 

t-T 

^(t-x) ^ 


w(x) dx (4.3) 


Since w(x) is a white noise, past measurements convey no infor- 
mation about it. Therefore j if no measurement is made between 
t^ and t, given an optimum estimate [x(t^)], the best estimate 
one can make about [x(t)] is 

[i(tn = [Jct*)] C4.4) 

(t. T>t>t.) 

^ 1 - 1 1 


a 


14 


The natural "damping" o£ the system tends to decrease the 
uncertainty in this estimate over that of the past estimate 
[x(tj^)]. However, the unknown white-noise excitation over 
this time period degrades the confidence one has in this 
new estimate. These two effects are uncorrelated and ad- 
ditive, resulting in a covariance matrix of the new estimate 
error given by 


[pct)]=Ecaict)]-[x(t)])c[;(t)]-[x(tD])''') 

“[♦(t.t.)] [P(t*)] (4.5) 

where the elements of [c()(t,t^)] and [Q(t,t^)] are given by 
Equations (3.7) and (3.13) with the quantity t there replaced 
by t-t^. 
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4.2 


IMPROVEMENT IN ESTIMATE RESULTING FROM A MEASUREMENT 


Consider an a priori unbiased estimate [x(tj)], 


together with its error covariance [P (t . ) ] =E ( [x(t . ) -x (t . ) ] 


J 


J 


[x(t j ) -x(t j ) ] ) are available prior to the next measurement 
An altimeter measurement 


y. = x^Ct.)+n. 


J 


(4.63 


E(n^) = 0 , E(np = R. 


is now made. Consider first, the altimeter is operating at 

low sampling rates. The measurement noise n may then be 

considered uncorrelated with past measurement noises and 

therefore uncorrelated with the a priori estimate error. 

One obtains immediately an improved a posteriori estimate 

of the geoid undulation x,(t.)> based on the measurement y. 

J - J 

and the a priori estimate x,(t.), as 

J 


^33 


^■^(t-) = x,(t.) + ^ (y.-x^(t-)) 


(4.7) 


Equation (4.7) follows from an application of Equation (4.0) 
and may also be considered the result of optimally combining 
two pieces of information (a priori and measurement) with 
uncorrelated errors. One may also rewrite this equation 
as , 


(S/N) 


X, (t - ) “X, (t') = "estimated signal" = (y.-x,(t.)) 

3 13 3 (S/N)+l J ^ 


(S/N) 


(S/N)+l 


("received signal") 
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A - A 

where (S/N) = P-^(t-)/R. = a priori signal- to-noise power 
ratio. Knowledge of an improved x,(t.) also enables us to 
obtain improved estimates of Xj^(tj) and X 2 (tj) based on the 
a priori correlation Pl 3 (tj) and PzsC^j)* Another application 
of Equation (4.0) with x^Ctp-x^CtT) playing the role of known 
quantity gives the a posteriori estimates 


x^Ctp-Xj^Ctp = 




(y.-XjCtJ) 


(4.8) 

4 


X2(t.)-X2(t.) = 4- 


P33(tp*Rj 




(4.9) 


Equations (4.7), (4.8) and (4.9) for the a posteriori estimate 
may be rewritten in the following compact matrix form 


[i(tp] = [i(t:)l + [K(tp] (y-J3(tp) 


(4.10) 


[K(t.)l = 


P33(t:).R. 




= optimum gain (4.11) 


The corresponding estimate errors may be written as 


[i(tp-x(t.)] = [i(t:)-x(tp] + [KCtpi (n-CijCtp-XjCt.))} 
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From this expression and noting that present measurement 
errors are uncorrelated with past estimates, one obtains the 
a posteriori estimate error covariance matrix 




= [PCtpi - [KCtp] [Pj3(t-):P23(t:);P33(t:)] C4.12) 


The optimum zero-lag filter may now be summarized as 
follows : 

1. Equations (4.4) and (4.10) are the recursive filter 
equations. 

2. Equations (4.5) and (4.12) are the recursive equation 
for the estimate error covariance. 

3. Equations (4.1) and (4.2) are initial conditions for 
starting these recursive equations. 

4. Equations (4.10) and (4.12) govern the interpolation 
between and extrapolation beyond measurements. 

If interpolation and extrapolation are not of 
interest, one uses in these equations to 

obtain a priori estimates just prior to the next 
measurements. 

5. In Equation (4.10) the difference between the new 
measurement and the a priori estimate is multiplied 
by a gain to obtain an optimum correction to the 


i 


a priori estimate. The optimum gain is computed 
from Equation (4.11), which shows its dependence 
on the a priori error covariance and the measure 
ment noise level. Equation (4.12) describes the 
decrease in the estimate error covariance as a 
result of measurement. 
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SECTION 5.0 
FIXED- INTERVAL FILTER 

5.1 OPTIMUM FIXED- INTERVAL FILTER* (SMOOTHER) 

The optimum fixed interval filter. Equation (2.5), 

makes use of all the altimeter data y^, i=0,l, . . .N-1. A 

zero- lag filter, on the other hand, uses only part of the 

data y., j=0,l,...q, (t £t=filter output time). However, 

once the optimum recursive zero-lag filter is constructed, 

an optimum smoother realizing Equation (2.5) may be obtained 

by properly combining a forward and a backward pass of the 

data through the zero-lag filter. Although this forward- 

backward-filter formulation of the optimum smoother is 

relatively recent (Reference 5) , the underlying idea is 

deceptively simple. The altimeter data yo’^l’*'’%-l 

first fed through the zero- lag optimal filter described in 

Section 4. The output of the filter gives us an estimate 

[x(t)] together with the error covariance [p(t)]. Since 

the geoid undulation is considered a stationary random 

process, which terminal of the data pass is the beginning 

and which in the end really does not matter. Therefore, 

one could feed the altimeter data in reverse order; i.e., 

i* 

^N-l’^N 2’'‘‘^0 same filter and obtain an 

estimate [x, (x)] and error covariance [Pv,(r)]. At any time 

D A, ^ D .A. 

t=t there exists two estimates, [x(t )] and [x- (x “Tj._^-t )] 
obtained by optimally processing altimeter data y^, i=0,l,...q 
and yj , j=N-l, N-2,...q+l, respectively. Since these esti- 
mates are based on different data which are corrupted only 
by unpredictable white noises, the estimate errors are 
uncorrelated. Therefore the over-all optimum estimate 
based on all the data may be obtained by optimally combining 
these two uncorrelated estimates. 


"Fixed-interval” refers to the fact that estimates of interest 
are to be derived from a common data set spanning a time interval 
of fixed duration. For conciseness, the terminology "smoother" 
will be used interchangeably. 

For non- stationary processes, a minor modification of the filter, 
reflecting the backward propagation of dynamics, is needed. 
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A word of caution is needed in combining the forward and back- 
ward filtered estimates because the forward and the backward 

state variables have different physical meanings. Although 
/\ /\ 

and x^^ represent estimates of the geoid undulations h, 
it is evident from Eq. (3.5) defining the state variables that 
'^2-a, estimate of (^+3)h while X2|^ is an estimate of (-^+3)h. 


if t is reckoned as the forward time 
to show that the transformation 


It is straight forward 


[x^] = [H][x^] 


[Pbl = [H][p^][H] 


with 


[H] 



-43 43^ 
-1 23 

0 1 


gives the backward estimates [x^^] and covariance [p^] which re- 
present the same quantities as the forward filter output. 

From the result given in the Appendix one obtains immediately 
the smoothed estimate by optimally combining the forward and 
backward estimates as 
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* - ” ^ * 


[x(t)] = [pCt)]_ ([p(t )] [x(t )] + [p, (t )] [x, (t )] , 




with an error covariance 


[P(t)] 


([pCt"" 


)] + [Pu(t-)] 


(5.2)’ 


Other theoretically equivalent but computationally more 

efficient ’algorithms are given in the Appendix. Notice' the 

smoother is symmetrical with respect to forward and backward 

filtering. It is immaterial whether forward or backward 

filtering is done first. If backward filtering is done first, 

'' + + 

the filter output )] and [P|^(t )] will be obtained to 

be stored in reverse order (normal order in t) . The quantities 
[x(t )] and [p(t ) ] are then computed "on-line” as altimeter data 
in normal order are fed-through, and these are combined "on- 

/s 

line" with [x^(t)] and [p^(t)] retrieved from storage to 
obtain the optimum smoothed estimate according to Equations 
(5.1) and (5.2). The major draw-back of this smoother is the 
complexity of data handling required; i.e., 

1. Reversing the ordering of altimeter data. 

2. Reversing the order of output of backward filter 


[^b^’^'i^^ and [P|^(t^)] and storing them for ready 


(3x1) 


(3x3, symmetric) 


retrieval later. 


In these formulas, the forward filtered estimate is combined 
with the backward one-step prediction in order to obtain the 
smoothed estimate. One may, of course, reverse the roles and 
combine the forward prediction with the backward filtering. 


I 


F 
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If one can tolerate the appearance of smoothed geoid undula- 
tions in reverse order, an alternative recursive smoother 
(Reference 6) may be used which operates on the outputs of 
the forward filter in a backward recursion. This eliminates 
the need for step ”1'* above, although operations equivalent 
to step ''2" above is still required. 

5.2 A PRIORI ACCURACY ASSESSMENTS AND A SUB-OPTIMAL SMOOTHER 

The covariance matrix gives us some information about 
the accuracy of the geoid undulation estimate. The equations 
governing the evolution of the covariance matrix in the pre- 
ceding sections show that the covariance matrix depends on our 
model of the signal and noise process, on measurement frequencies, 
etc., but does not depend on the actual altimeter data. This 
means that if our model is reasonably accurate, the accuracy of 
the geoid undulation estimate may be assessed beforehand by 
solving the variance equations. The results form a standard 
by which one may measure the performance of simpler suboptimal 
smoothers. Also of interest are the convergence characteristics 
of the variance equation. 

The optimum smoother. Equations (5.1) and (5.2), requires 

the input of the filtered state and covariance obtained during 

the first filter pass. In addition to the storage of the 

filtered outputs for retrieval in reverse order, inversions 

of the covariance matrices at each measurement time are 

required. Since the dimension of the geoid undulation state 

is 3, this portion of the smoother alone requires the reverse 

storage of (3 state variables + 6 independent elements of 

covariance matrix) * (N measurements) and the inversion of a 

t 

3x3 covariance matrix N times. In order to reduce the com- 
putational and storage requirements, the following sub-optimal 

•*. 

Equations (5.1) and (5.2) as they stand require 3N inversions. 

Other algorithms given in the Appendix show this may be 
reduced to N inversions. 
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smoother may be suggested. The forward and backward filtering 
per se do not involve particular difficulties and will be 
implemented. The optimal processing of the two filtered 
estimates to obtain the optimum smoothed estimate as de- 
scribed by Equations (5.1) and (5.2) will be replaced by their 
scalar versions involving the geoid undulation only; i.e., 
replacing [x] by and [p] by in these equations. It is 

expected that this sub-optimal estimate will be reasonably 
good for the following reason. Although there is a degrada- 

tion in the estimate (x^(t)) resulting from not utilizing 

^ s 

information about the other state variables, yet this degrada- 

ys 

tion does not propagate in time because (x,(a)) for a^t does 

y\ O 5 

not depend on (x,(t)) . In any case, the following procedure 

s 

may be suggested for assessing the quality of this sub-optimal 
smoother. The estimate covariance at a few selected instants 
during the data pass may be computed from Equation (5.2). 

Since it is only for a few points, the computational and 
storage requirements are minimal. These minimum geoid 
undulation variances may be compared with those obtained 
from the sub-optimal filter. In particular, one should 
expect the best accuracy exists near the middle of the data 
pass and the minimum variance obtained there provides a lower 
bound for the variance, or, an indication of the achievable 
accuracy. 

To illustrate the ideas discussed above, let us consider 

a typical GEOS- III pass with altimeter measurements sampled at 

a nominal rate of 10 measurements per second (sampling interval 

= 0.102406 sec.). The geoid undulation is considered to have 

2 

an a priori variance of 4 meters and a correlation length of 
50 Km. (g=0 . 380 5/sec. ) . The altimeter noise is assumed uncor- 
related from sample to sample and to have a variance of 0.36 

2 • 
meter , corresponding to an input signal to noise power ratio 

of 4/0.36 « 10. Solutions to the variance equations give the 

following results: 
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Forward filtered error covariance matrix reaches 
a steady-state value of 

0,1528646 

0.9721489x10"^ 0,1036214 

0.4034996x10"^ 0.6118031x10"^ 0.5292089x10" 

Backward filtered one-step prediction error covari 
ance matrix (untransformed) reaches a steady-state 
value of 

0.1581666 

0.1052539 0.1158105 

0.4730372x10"^ 0.7172390x10"^ 0.6204108x10" 

Forward-backward smoothed error covariance matrix 
reaches a steady-state value of 

0.3134349x10"^ 

0.6365069x10"^ 0.1338882x10"^ 

-0.6004677x10"^ 0.5315084x10"^ 0.1396868x10" 

One-way forward filter gain reaches a steady- state 
value of 

“0.1120832" 

0.1699453 

0.1470024 

In the steady-state, the contributions of the 
forward filtered and backward one- step predicted 
estimates to the smoothed estimates may be seen 
by the respective weighting matrices, i.e., 


Forward 


" 0.3486159 

-0.1313888 

-0.2273753 

Weighting 

[W]p = 

-0.1989280 

0 . 5260603 

-0.3560544 

Matrix 


0.2653295x10"^ 

-0.3293662 

0.6447248 

^ r* Vw?? 


" 0.6513841 

-0.1122795x10 

0.7045924 

CL i\. W CX X \J. 

Weighting 

[W]b = 

0.1989280 

-0.7767080 

0.8319259 

Matrix 


-0.2653295x10"^ 

-0.3293661 

0.6059229 


Note [W]g is the weighting matrix for the un- 
transformed backward estimate; i.e., [x^] in 

Section 5.1. The weighting matrix for the 
transformed backward estimate [x^] is, of 
course, the difference between the identity 
matrix and [W]p. It may be seen from these 
weighting matrices that the direct contributions 
of the estimates of x^ to the smoothed value of 
the geoid undulation are insignificant. 

6. The convergence characteristics of the one-way 
filter are described by the graphs in Figure 5.1. 

It is seen that the filter pretty much reaches 

the steady state after 30 measurements (3 seconds) . 

7. For the sub-optimal smoother described previously, 
the standard deviation of the geoid undulation 
reaches a steady-state value of /0 . 2 85 5961x10"^ 

= 0.169 meter. This compares with the optimal 
value of /O . 1396868x10"! = 0.118 meter for the 
optimal smoother. 

8. Sometimes, in addition to geoid undulation, the 
geoid undulation slope is also of interest. The 
slope = Xp-8 Xj, and simple calculation shows 
the smoothed undulation and slope covariance 
matrix reaches a steady-state value of 


V 


I 


0.1136643x10"^ 0(10"^) 

0(10'^) 0.1396869x10 

It is o£ interest to note the smoothed geoid 
undulation and slope are almost uncorrelated. 

The slope estimation accurac y is specified 
by a standard deviation of *^0.1136643x10"^ 

= 0.1066135 meter/sec. =: 3.36 arc sec. at the 
approximate GEOS-III ground speed of 6.55 Km/sec, 
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GO UP TO A PRIORI VALUE 
OF 2 METERS 


FIGURE 5.1. 

CONVERGENCE CHARACTERISTICS OF ONE-WAY 
KALMAN FILTER 
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5.3 ASYMPTOTIC BEHAVIOR OF SMOOTHER- FREQUENCY RESPONSE 

AND DATA WEIGHTING SEQUENCE 

Away from the beginning and the end of a long data 
sequence the smoother is governed essentially by its steady- 
state asymptotic behavior. The discrete transfer function 
describing the smoother asymptotic characteristics is derived 
in Appendix B. From the transfer function, quantities which 
give one intuitive understanding of the filter such as the 
frequency response and the data weighting sequence may be 
calculated immediately. To illustrate, the same numerical 
example as used in Section 5.2 will be considered. Substi- 
tution of the appropriate steady-state smoother parameters 
given in Section 5.2 into Equation (B.5), and after some 
algebra, one obtains the following transfer function relating 

/s 

x^, the geoid height estimate, to the altimeter measurement y. 


H(z) = 0.0388018553 + f(z) + f (i) 


where 

0 .0386814002 z^- 0,06 71610539 z+0. 029446601 7 

f(z) = -7 7 

z^-2. 7 26669377z^+2. 487576809 z-0. 7588967469 


The frequency response of the smoother is obtained 
from the transfer function by replacing z by i.e., 

HCe^“'^) = 0.0388018553 + fCe’’^'^) + f(e~^'^'^) 

= 0.0388018553 + 2 Re {f(e^^'^)} 
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That the above expression is real indicates there is no phase 
lag, a fact which could have been anticipated from the com- 
pletely symmetrical forward and backward operation of the 
smoother. The amplitude response as a function of excitation 
frequency is computed, tabulated in Table 5.1 and plotted in 
Figure 5.2. 

The data weighting sequence, or the variations of data 
weights with data "distance" from the estimate are coefficients 
of the polynomial expansion of the transfer function. As ex- 
pected, because the dependence of H(z) on "z" and "1/z" is 
identical, there is no preference to "past" or "future" data. 

The data weighting sequence is computed, and is tabulated in 
Table 5.2, and plotted in Figure 5.3. Note that data weights 
exhibit a damped oscillatory behavior with a period of approxi- 
mately 57 sampling intervals (5.8 seconds, or 38 Km). It is 
also of interest to note that as shown in the preceding section, 
the current measurement is given the weight 0.147 by the forward 
filter. With the smoother, the weight assigned is 0.0388, 
indicating the much improved smoothing effect. 
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FIGURE 5.2. SMOOTHER FREQUENCY RESPONSE 


i 



A Priori Signal to Noise Variance Ratio = 11.1 
Signal Correlation Length = 50 Km 


• In Phase 
o 180° Out of Phase 




Table 5.1. Smoother Frequency Response 


FREQUENCY, 

0 

0,0326 

0.0651 

0.0976 

0.130 

0.163 

0.195 

0.288 

0.260 

0.293 

0.326 

0.358 

0.391 

,0.423 

0.455 

0.488 


0.520 


-74.0 

0.553 


-66.3 

0.586 


-64.4 

0.618 


-63.9 

0.651 


-64.1 

0.683 


-64.6 

1.01 


-72.7 

2.02 


-95.6 

3.03 


-109. 

4.88 

(Nyquist Frequency) 

-120 


AMPLIFICATION, db 


+ 0.005 
-0.037 
-0.190 
-0.580 
-1.54 
-3.64 
-7.29 
-12.3 
-18.0 
-24.0 
-30.0 
-36.0 
-42.1 
-48.7 
-56.4 
-68.4 


There is a 180° phase shift for responses to frequencies 
higher than =**0.5 Hz 
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DATA WEIGHT 
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figure 5.3. DATA WEIGHT VS DATA DISTANCE FROM ESTIMATE 



(SAMPLING PERIOD = 0.1024 *ec., equivalent to 0.67 Km. ground distance) 


Table 5.2. Variation of Data Weights With 
"Distance" From Estimate 
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DATA "DISTANCE" FROM ESTIMATE, 

DATA WEIGHTS 

IN TERMS OF NUMBER OF SAMPLING 
PERIODS, N 

(Sampling Period=0 . 1024 sec. 
equivalent to 0.67 Km of 
ground distance) 

INSTANTANEOUS 

W(N) 

CUMULATIVE 

N 

W(0)+2 ^ W(n) 
n=l 

0 

0.388019D-01 

0.038802 

1 

0.386814D-01 

0.116165 

2 

0. 383103D-01 

0.192785 

3 

0.376833D-01 

0.268151 

4 

0.368051D-01 

0.341762 

5 

0. 356889D-01 

0.413140 

6 

0.343540D-01 

0.481848 

7 

0.328245D-01 

0.547497 

8 

0.311274D-01 

0.609752 

9 

0.292918D-01 

0.668335 

10 

0. 273478D-01 

0.723031 

11 

0.253252D-01 

0.773681 

12 

0. 232531D-01 

0.820188 

13 

0.211594D-01 

0.862506 

14 

O.19O7O0D-O1 

0.900646 

15 

0.170086D-01 

0.984663 

16 

0.149965D-01 

0.964657 

17 

0.130526D-01 

0.990762 

18 

0.111929D-01 

1.013147 

19 

0.943070D-02 

1.032009 

20 

0.777686D-02 

1.047563 

21 

0.623958D-02 

1.060042 

22 

0.482466D-02 

1.069691 

23 

0.353564D-02 

1.076762 

24 

0.237402D-02 

1.081510 

25 

0. 133940D-02 

1.084189 

26 

0.429734D-03 

1.085049 

27 

-0. 358482D-03 

1.084332 

28 

-0 .102999D-02 

1.082272 

29 

-0 . 159057D-02 

1.079091 

30 

-0. 204683D-02 

1.074997 

31 

-0 . 240602D-02 

1.070185 

32 

-0.267585D-02 

1.064833 

33 

-0.286433D-02 

1.059104 

34 

-0.297963D-02 

1.053145 

35 

-0 . 302990D-02 

1.047085 

36 

-0 . 302323D-02 

1.041039 

37 

-0.296745D-02 

1.035104 

38 

-0 . 287013D-02 

1.029364 
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Table 5 . 2 (Cont * ) 


DATA "DISTANCE" FROM ESTIMATE, 
IN TERMS OF NUMBER OF SAMPLING 
PERIODS, N 

(Sampling Period=0 . 1024 sec. 
equivalent to 0.67 Km of 
ground distance) 

DATA WEIGHTS 

INSTANTANEOUS 

W(N) 

CUMULATIVE 

N 

W(0) + 2 ^ W(n) 
n=l 

39 

-0.273846D-02 

1.023887 

40 

-0.257919D-02 

1.018729 

41 

-0.239859D-02 

1.013931 

42 

-0.220246D-02 

1.009526 

43 

-0.199602D-02 

1.005534 

44 

-0-178399D-02 

1.001966 

45 

-0.157054D-02 

0.998825 

46 

-0-135930D-02 

0.996107 

47 

-0.115339D-02 

0.993800 

48 

-0.955425D-03 

0.991889 

49 

-0.767552D-03 

0.990354 

50 

-0.591470D-03 

0.989171 

51 

-0.428469D-03 

0.988314 

52 

-0.279458D-03 

0.987755 

53 

-0.145005D-03 

0.987465 

54 

-0.253708D-04 

0.987414 

55 

0.794530D-04 

0.987573 

56 ^ 

0.169710D-03 

0.987913 

57 

0.245845D-03 

0.988404 

58 

0.308466D-03 

0.989021 

59 

O.358321D-03 

0.989738 

60 

0.396259D-03 

0.990531 

61 

0.423211D-03 

0.991377 

62 

0.440161D-03 

0.992257 

63 

0.448123D-03 

0.993154 

64 

0.448121D-03 

0.994050 

65 

0.441176D-03 

0 .994932 

66 

0.428284D-03 

0.995789 

67 

0.410408D-03 

0.996610 

68 

0.388463D-03 

0.997386 

69 

0.363314D-03 

0.998113 

70 

0.335762D-03 

0.998785 

71 

0.306543D-03 

0.999398 

72 

0.276327D-03 

0.999950 

73 

0. 245710D-03 

1.000442 

74 

0.215222D-03 

1.000872 

75 

0 .185318D-03 

1.001243 

76 

0.156390D-03 

1.001556 

77 

0 .128762D-03 

1.001813 

78 

0 .102696D-03 

1.002019 

79 

0.783964D-04 

1.002175 

80 

0 560145D-04 

1.002287 
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Table 5.2 (Cont.) 


DATA "DISTANCE" FROM ESTIMATE, 

DATA WEIGHTS 

IN TERMS OF NUMBER OF SAMPLING 
PERIODS, N 

(Sampling Period=0 . 1024 sec. 
equivalent to 0.67 Km of 
ground distance) 

INSTANTANEOUS 

W(N} 

CUMULATIVE 

N 

W(0) + 2 ^ W(n) 
n=l 

81 

0.356515D-04 

1.002359 

82 

0.173640D-04 

1.002393 

83 

0.116954D-05 

1.002396 

84 

-0.129497D-04 

1.002370 

85 

-0.250413D-04 

1.002320 

86 

-0.351785D-04 

1.002249 

87 

-0.434554D-04 

1.002162 

88 

-0.499830D-04 

1.002063 

89 

-0. 548854D-04 

1.001953 

90 

-0.582960D-04 

1.001836 

91 

-0.603541D-04 

1.001715 

92 

-0.612023D-04 

1.001593 

93 

-0.609836D-04 

1.001471 

94 

-0.598392D-04 

1.001351 

95 

"-0.579066D-04 

1.001236 

96 

-0*.553178D-04 

1.001125 

97 

-0.521980D-04 

1.001021 

98 

-0.486645D-04 

1.000923 

99 

-0.448261D-04 

1.000834 

100 

-0.407820D-04 

1.000752 

101 

-0.366221D-04 

1.000679 

102 

-0.324263D-04 

1.000614 

103 

-0.282649D-04 

1.000557 

104 

-0.241985D-04 

1.000509 

105 

-0.202784D-04 

1.000468 

106 

-0.165471D-04 

1.000435 

107 

-0.130384D-04 

1.000409 

108 

-0.977857D-05 

1.000390 

109 

-0.678637D-05 

1.000376 

110 

-0.407406D-05 

1.000368 

111 

-0.164793D-05 

1.000365 

112 

0.491027D-06 

1.000366 

113 

0. 234643D-05 

1.000370 

114 

0. 392586D-05 

1.000378 

115 

0.524024D-05 

1.000389 

116 

0.630322D-05 

1.000401 

117 

0.713061D-05 

1.000416 . 

118 

0.773989D-05 

1.000431 

119 

0. 814967D-05 

1.000447 

120 

0. 837927D-05 

1.000464 

121 

0. 844836D-05 

1.000481 

122 

0. 837656D-05 

1.000498 

123 

0. 818316D-05 

1.000514 

124 

0. 788688D-05 

1.000530 

125 

0.750561D-05 

1.000545 
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Table 5.2 (Cent.) 


DATA "DISTANCE” FROM ESTIMATE, 

DATA WEIGHTS 

IN TERMS OF NUMBER OF SAMPLING 
PERIODS, N 




CUMULATIVE 

(Sampling Period=0 . 1024 sec. 

INSTANTANEOUS 

N 

equivalent to 0.67 Km of 
ground distance) 

W(N) 

W(0) + 2 \ W(n) 


n=l 

126 

0 .705627D-05 

1.000559 

127 

0.655467D-05 

1.000572 

128 

0.601538D-05 

1.000584 

129 

0.545169D-05 

1.000595 

130 

0.487555D-05 

1.000605 

131 

0.429758D-05 

1.000613 

132 

0.372703D-05 

1.000621 

133 

0.317186D-05 

1.000627 

134 

0.263877D-05 

1.000632 

135 

0.213323D-05 

1.000637 

136 

0.165959D-05 

1.000640 

137 

0.122112D-05 

1.000642 

138 

0.820151D-06 

1.000644 

139 

0.458100D-06 

1.000645 

140 

0.135603D-06 

1.000645 

141 

-0.147403D-06 

1.000645 

142 

-0.391592D-06 

1.000644 

143 

-0.598158D-06 

1.000643 

144 

-0.768726D-06 

1.000642 

145 

-0.905276D-06 

1.000640 

146 

-0.101006D-05 

1.000638 

147 

-0.108555D-05 

1.000636 

148 

-0.113434D-05 

1.000633 

149 

-0.115910D-05 

1.000631 

150 

-0.116257D-05 

1.000629 

151 

-0.114742D-05 

1.000626 

152 

-0.111629D-05 

1.000624 

153 

-0.107175D-05 

1.000622 

154 

-0.101620D-05 

1.000620 

155 

-0.951940D-06 

1.000618 

156 

-0.881099D-06 

1.000616 

157 

-0. 805631D-06 

1.000615 

158 

-0. 727312D-06 

1.000613 

159 

-0.647733D-06 

1.000612 

160 

-0.568301D-06 

1.000611 

161 

-0.490238D-06 

1.000610 

162 

-0.414586D-06 

1.000609 

163 

-0. 342217D-06 

1.000608 

164 

-0. 273838D-06 

1.000608 

165 

-0. 210002D-06 

1.000607 

166 

-0.151121D-06 

1.000607 

167 

-0.974754D-07 

1.000607 

168 

-0.492279D-07 

1.000607 

169 

-0.643593D-08 

1.000607 

170 

0. 309357D-07 

1.000607 

171 

0.630025D-07 

1.000607 

172 

0.899478D-07 

1.000607 

173 

0. 112011D-06 

1.000607 
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/ 




DATA "DISTANCE" FROM ESTIMATE, 

DATA WEIGHTS 

IN TERMS OF NUMBER OF SAMPLING 
PERIODS, N 


CUMULATIVE 

(Sampling Period=0 . 1024 sec. 

INSTANTANEOUS 

N 

equivalent to 0.67 Km of 

W(N) 

W(0)+2 y ^ W(n) 

ground distance) 


n=l 

174 

0.129478D-06 

1.000608 

175 

0.142669D-06 

1.000608 

176 

0.151928D-06 

1.000608 

177 

0.157620D-06 

1.000608 

178 

0.160114D-06 

1.000609 

179 

0.159785D-06 

1.000609 

180 

0.157001D-06 

1.000609 

181 

0.152124D-06 

1.000610 

182 

0.145499D-06 

1.000610 

183 

0. 137456D-06 

1.000610 

184 

0.128303D-06 

1.000610 

185 

0.118326D-06 

1.000611 

186 

0.107787D-06 

1.000611 

187 

0.969242D-07 

1.000611 

188 

0.8 59 48 3D- 0 7 

1.000611 

189 

0.750457D-07 

1.000611 

190 

0.643771D-07 

1.000612 

191 

0. 540791D-07 

1.000612 

192 

0. 442648D-07 

1.000612 

193 

0.350250D-07 

1.000612 

194 

0. 264301D-07 

1.000612 

195 

0.185312D-07 

1.000612 

196 

0.113618D-07 

1.000612 

197 

0.493985D-08 

1.000612 

198 

-0.730738D-09 

1.000612 

199 

-0. 565831D-08 

1.000612 

200 

i 

-0.986174D-08 

1.000612 







5.4 EFFECT OF INCORRECT GEOID UNDULATION MODEL 

The geoid undulation has been considered as a sta- 
tionary random process with a zero mean value and with a 
TASC correlation function. Equation (3.1). The TASC model 
is but one of many possible correlation functions which 
satisfy certain very general criteria for geoid undulations 
(Ref. 2) . There is no a priori reason to believe the TASC 
correlation function is the model, much preferable to others. 
Instead one may adopt the point of view that the TASC model 
is one of many reasonable models based on which a filter 
may be designed rationally with the desirable features of 
being simple, symmetric and having fading memory. The two 
parameters which define the TASC correlation function are: 

2 

1. The variance a . 

2. The correlation length, or equivalently the 
"time constant" 4, 

p 

The variance equation and the optimum gain given in Section 4 

may be normalized with respect to the geoid undulation vari- 

2 2 
ance a . It may be seen then that instead of a , it is the 

2 

signal- to-noise ratio a which is of importance. The only 
direct effect of S/N ratio is on the gain. Naturally the 
desired gain increases and decreases with the increase and 
the decrease of the signal- to-noise ratio. The correlation 
length, or "time constant" 4 , determines the effective 

p 

memory length of the filter, and if it is taken to be too 
large, sharp features will be smeared. On the other hand, 
if the correlation length is too short, the filter cannot 
differentiate the higher frequency noise from the lower 
frequency geoid undulation. Although the effects of the 
two parameters are somewhat coupled, broadly speaking, when 
examining the smoothed altimeter data, if one feels the 
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amplitudes are too low, the S/N ratio assumed is probably 
too small; if the features are not sharp enough, the cor- 
relation length is likely too long. 

Both the signal and noise are assumed to be zero 
mean random processes. The actual data may contain low 
frequency trends due to biases, orbital errors, and etc. 
Ideally, if geoid undulation is the output desired, these 
trends should be removed before the data are input to the 
filter. Sometimes, the operational objective is just to 
"smooth" the data, leaving the trends intact. Generally 
speaking, filters should have "low pass" features and 
we shall investigate the ability of our smoother in passing 
through the low frequency trends. The signal and noise 
models shall be the same as those considered in the pre- 
ceding section. The trends considered are a constant bias 
and a constant slope. The following table summarizes the 
steady-state* smoother outputs resulting from these trends. 
The results, of course, agree with the frequency response 
given in the preceding section. 


Away from the beginning and the end of the data pass . 
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SIMULATION RESULTS 


TREND 

FORWARD 

ESTIMATE 

BACKWARD I 

PREDICTION 

SMOOTHED VALUE 


HEIGHT 

m. 

SLOPE 
m/sec . 

HEIGHT 

m. 

SLOPE 
m/sec . 

HEIGHT 

m. 

SLOPE 
m/sec . 

+ 1 Meter 
BIAS 

0.976325 

-0.038262 

0.972245 

0.039874 

1.000611 

0.000000 

ZERO HEIGHT 

0.097748 m/sec, 
SLOPE 

-0.014554 

0.073774 

0.012844 

0.072628 

0.000000 

0.091840 


























1 


J 


I 


I 


w 


i 


Ik 


The simulation results show that a constant trend or 
bias in the data will have virtually no effect on the slope 
estimate but will cause a 0.06% height error. On the other 
hand, a linear trend, or constant drift in the data will have 
no effect on the height estimate but will cause a 6% slope 
error . 


Note that the forward filter ''underestimates” the height 
and slope of the bias. This is not unexpected since the signal 
is modeled as a stable system,* and will be damped between 
measurements. In particular, it will be damped to zero for 
a large data gap. This is not unreasonable as it says the 
geoid undulation estimate should be zero since the geoid 
undulation is just as likely to be positive as negative in 
the absence of any measurement information. That both the 
height and the slope are under-estimated is consistent with 
the positive correlation between these quantities as indicated 
by the forward filtered error covariance matrix given in 
Section 5.2. The backward filter exhibits a similar behavior 
to the bias. It underestimates the height and the "backward 
slope." But since the backv;ard slope is the negative of the 
forward slope, the backward filter over-estimates the forward 
slope and also exhibits a negative correlation between the 
height and the "forward slope estimates." It is of interest 
to note that although both the forward and backward filters 
under-estimate the height, the differences in their slope 
estimates compensate each other and the resulting optimum 
smoothed results not only give an almost perfect slope esti- 
mate, but also bring along with it a slightly over-estimated 
height. A somewhat analogous situation exists for the linear 
trend. 


Tracking unstable systems with accuracy is generally impracti- 
cal. For a discussion of that subject see Ref. 8. 
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SECTION 6.0 

OPTIMUM RECURSIVE FILTER AT HIGH ALTIMETER SAMPLING 
RATES -MODIFICATIONS FOR CORRELATED NOISE SEQUENCE 


Some correlations may exist among the altimeter noises 
when measurements are made at rates greater than 10 samples/ 
sec. Modifications of the results of Section 4.2 to take 
into account this correlation are presented in this section. 
For a general treatment of correlated noise see Reference 7. 

Consider altimeter measurements at two instants, t- 

’ 1 

and tj ; i . e . , 


= XjCt.J + n. 

(6.1) 

yj = 

(6.2) 


The existence of correlation between the noises n- and n- 
means knowledge of n^ conveys information about n^ , or, 
past measurements enable us to improve our a priori knowledge 
of future noises. As a matter of fact, an application of 
Equation (4.0) enables us to relate n^ and n. in the follow- 
ing way, 


^ ^ A E(n.-n.) 

n. = n. + (n.-n.) = ■■ ■ ^ n. + u 

J ^ J ^ E(n.) 




(6.3) 


Notice that Equation (3.17), given before, is just a special 
case of Equation (6.3) for exponentially correlated noise, 
which is our model of the altimeter noise at high sampliiVg 
rates . 


By making use of Equations (3.17) and (6.1) and taking 
j=i+l, one may replace Equation (6.2) by 


y 


bCti.i-ti) 


i+1 






X3(t.)*e 


b(ti^l-tp 


u . 

1 


Furthermore, because geoid undulations at "adjacent locations" 
are related by Equation (3.8), one obtains finally 




1 - 



[x(t^] 

3x1 



-B(Vi-x) 


w (t) 


bTi 

dx+e u^. 


where = ^i+l~^i sampling period. The above equation 

is in the standard form of the measurement equation in the 
Kalman filter; i.e.. 


■i = [H(t.)] [x(t.)J + V. , 


C6.4) 


with 



bT- 



e ^ 


_ 

BT. T 

2 

[H(t^)] = e 

^ [- 

f ’’i 1 

V. 

1 A. 
1 

(*■ - t) ^ 

^"i+1 ^ 

-ecti.i 

A 

2 


E(v.) = 0 , 

E(v.v.) 

= G. 6 . . 

1 ij 

2bT. 


Gi = Ri te 

^-1) + q 

33 


BT. 


I bT. 

w(t) dx+e ^ u> 
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With Equation (6 . 4) replacing the original measurement equation, 
a Kalman Filter for the correlated altimeter data may be con- 
structed immediately. Note the following complexities are 
introduced because of the altimeter noise correlations: 

1. The original altimeter data have to be converted 
to an equivalent data set. This involves some 
computations and introduces a one-sample-period 
lag to the filter output. 


2. The converted measurements are related to all 

three state variables instead of just the geoid 
undulation. 


3. Correlations now exist between the measurement 


noise and the dynamic noise; i.e., 


•ti+i -3(t. -,-t) 
Efv. / "-^^e 


/ 




w(t) dr 


t- 


,2 r 

t. 

1 


^13 
*^23 
533 J 






dr 


Although "2” and "3” do not present any difficulties, they 
do make the filter structure more complex. 
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SECTION 7.0 

PROGRAM IMPLEMENTATION 


A computer program, ALTKAL, which inqjianents the algoriths 
described in the preceding sections has been developed at the Wallops 
Flight Center, Wallops Island, Virginia. The program is designed 
specifically for the smoothing of sea surface heights derived from GEOS-3 
altimeter data. The program ccm^utes, and outputs, filtered sea surface 
heights from both a forward and a backward pass through a data set. The 
forward and backward filter results are then ccmibined to obtain the optimally 
smoothed heights. The con 5 >utations include the smoothed sea surface slope, 
which is also output by the program. 

Although ALTKAL was intended to be used primarily with cumulative 
altitude data (data period of 0.1024 seconds), any regular or irregular data 
rate can be used. However, the iiqjlemented algorithms assume that measure- 
ment noise is uncorrelated between san^les, so correlated input (such as 
instantaneous altitudes from GEOS-3) would not be optimally treated. 

7.1 PROGRAM INPUT 

The program will accept as input either of two types of input 
tapes. These tapes are: 

GEODYN binary residual tape 

ARC tape 

Detailed formats for these tapes are given in references 9 and 10, 
respectively. 
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7.2 PROGRAM OUTPUT 

The complete output o£ the program includes : 

• A list of the input parameters used in the 
TASC model. 

• A complete listing of the processing results. 

• The RMS of the residuals of the measurements 
about the final smoothed heights. 

• An output tape containing relevant variables 
to be saved or for use in plotting (optional) . 

• An output tape suitable for input to the 
SEAHT program (optional) . 

Specific examples and discussion of the output can be found 
in Section 8.3. 

7.3 OPTIONS 

ALTKAL has the capability of allov-zing specification 
of the input processing paYameters used. Through the- use 
of keyword option cards (Section 8.2) certain options can 
be implemented. These include; 

• Adjusting the editing threshold. 

• Adjusting the correlation length. 

• Adjusting the geoid undulation sigma. 

• Adjusting the altimeter noise sigma. 

• The manual editing (culling) of data points. 
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SECTION 8.0 
PROGRAM OPERATION 


( 


The operational ALTKAL program resides as a card deck 
at the Wolf Research and Development Group office in Pocomoke 
City, Maryland. The deck is operational on the Honeywell 
625/35 system at NASA/Wallops Flight Center. The run deck 
consists of the appropriate control cards, the program 
(binary cards), and the desired keyword cards. The standard 
set up for a GEOS-3 data run consists of the following: 

• Specifying the FORTRAN logical units used for 
input and output procedure. 

• Designating the appropriate options desired. 

A specific example of a complete set up can be found 
in Section 8.3.1. 


s 


8.1 FORTRAN LOGICAL UNITS 

FORTRAN 

LOGICAL 

UNIT USE 

1 GEODYN binary residual tape input^ 

2 ARC tape input^ 

11,12 Disk storage used during execution 

* Output tape used for plotting 

** Output tape used as input to SEAHT program 


The input tape number on which the individual pass is 
found must be specified on either 1 or 2 as indicated on 
the DATA option card. 


This can be any unit except 1, 2, 11, or 12 as specified 
on the PLOT option card. 


This can be any unit except 1, 2, 11, or 12 as specified 
on the SEAHT option card. 
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8.2 KEYWORD OPTION CARDS 

Each option card is identified by the name beginning 
in column 1 and is read under an A6 format. The name is 
followed by 5 fields: 

Field Columns Format 


1-6 

A6 

7-10 

14 

11-25 

E15.6 

26-40 

E15.6 

41-55 

E15.6 

56-70 

E15.6 


The option cards can be placed in any order and, if the card 
is omitted, all option variables have default values as 
indicated . 

The option cards are: 

CULL 

CORLEN 

DATA 

EDIT 

PLOT 

REVNUM 

SEAHT 

SIGMA 


50 


CORLEN 


CDRLENj 


000000 
1 } : < > (l 
1 1 1 1 n| 


2 2 2 2 2 

3 3 3 

4 4 4 4 4 

5 5 5 5 

6 6 0 0 

? 7 7 7 7 )| 


0 0 0 0 0 

I J II II 
11111 


2 2 2 2 2 2 
3 3 3 3 3 3 
44 4 4 4 4 
5 5 5 5 5 
6|B 6 6 6 6 
7 7 7 7 7 


25 . 


0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

17 17 II IS l( II II 19 70 71 77 73 71 7 S 71 71 71 79 II 31 II II II IS It II II I! ID ll II II II IS IS II II 19 SI SI 57 SI H 55 5 E 51 51 59 SO II SI II SI SS SI El E 3 tl ID II 17 73 II IS IE II II IS ID 

111111111111111111111111111111111111111111111111111111111111111111111 

22222222222222222222222222222222222722222222222222222222222222222272 

3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 ,1 

444444444444444444444444444444444444444444444444444444444444444444444 

5 5555555555555555555555555555555555555555555555555555555555555555555 

6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 E 6 O b 6 6 6 6 6 6 6 6 6 6 6 6 6 6 0 6 6 6 6 6 6 6 u 6 6 6 6 6 6 6 6 6 6 6 6 6 E 6 6 6 6 6 6 8 6 6 6 6 e E 

7 7 7 7 7 7 7 7 7 7 7 7 ’ 7 7 7 7 7 7 7 7 7 7 7 7 7 n 7 7 I 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 1 7 7 7 7 7 7 7 7 7 7 

8 8 8 8 8 8 I 8 8 8 8 8 8 8 8 8 8 8 3 8 8 0 e e 3 6 6 8 S B e 6 8 8 8 8 6 8 8 8 8 8 8 e 8 8 8 ! 0 8 S 8 8 8 I 8 8 8 1 8 8 J 8 B e 8 8 8 

3 9 9 9 JS 9 9 S|9 9 9 9 9 9 9 9 3 9 9 9 9 9 9 9 9 9 3 9 3 9 9 9 3 9 9 9 S 9 9 9 S 5 9 9 9 9 9 9 9 9 9 8 8 9 9 9 9 9 9 9 9 9 9 9 9 S 9 9 9 9 9 9 9 !> 9 9 9 3 

I 1 ? 1 } M U n II IJ ?» n !) H li K V 7 t 7 ! ir 37 J 3 V Ji 3 b lif If ‘‘ •? <1 H lb 1 j V U I 8 .SO b« b 7 M 51 55 55 5 » 'A 55 60 5 I 6 ? 6 ] 51 6 b 66 6 > 51 5 ! 70 7 t I? 71 H 15 75 JI H 75 IB 
MECC-bOBi 


NAME 

COLUMN 

FORMAT 

DESCRIPTION 

CORLEN 

1-6 

A6 

Modifies correlation length used 
in TASC model. 


11-25 

E15.6 

Correlation length in meters = 
CORLEN. 

IF CARD 

OMITTED: 

Default 

value of CORLEN = 50.0 is used. 


CULL 


CULL 


I ooooDOBOODOooDqcoosoitoenooBiiaiiioaoiiiitgooDoooooogDMooBOOooooieoiesiiiiPioii 

I i : i i i - I 1 It II li |] M IS K IMI I! il II II :i It IS » nil I) U II U U 34 U » )l U is 4t 41 41 II 14 IS II 41 41 I! Si 31 SI S3 i> SS SI SI S) S3 II II II II II IS H II II IS II II II II |4 IS II II II II U 

I 11 I II n 1 1 1 M t I I II M n I 11 11 11 11 I I M M 11 I II I I 1 1 1 1 M 11 1 1 I I 1 1 II I 1 1 1 I 1 1 11 I I I I 1 I 1 1 I 

2222222222222 222222222222222222222222222222222222222222222222222222222222222222 
3 3333333333 3333333333333 33333333333333333333333333333333333333333333333333 

4 444444444444444444444444444444444444444444444444444444444444444444444444444444 

5 S 5 3 5 5 S S S S 5 5 S S S 5 S 3 S S S 5 S S 5 5 5 5 5 5 5 5 5 5 S 5 S 5 5 S S 5 5 S S S 5 5 5 5 5 5 5 S 5 S 5 3 S 5 5 5 S S S 5 S S 5 S 5 5 5 3 S 5 S 5 5 

I i a $ I a S 6 t 6 E 6 6 6 S S S a 6 E 6 6 6 6 i 6 B S 0 b G S 6 S 6 6 6 6 6 ( 6 6 6 S S 6 6 E S & { 6 S 8 8 6 S B 6 a a S E fi t i i C S 6 ( a 5 G 6 6 S E S S 

1 n n m m n 111 m 1 1 11 :i n i n 1 1 1 1 n n 1 1 1 1 1 1 1 1 n 1 1 1 1 1 1 1 1 n 1 1 1 1 1 n 1 1 n i n 1 1 n 1 1 1 

a I a I a a a a I a 1 1 a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a 
saasss!S9sassss99iassa!ss!sssMaiaaasisaaaaaass99S98S99sa9999999saaas!sssasssan 

I I I I s I ' I ) <1 I' I II It IS II II II I] II ii II II II II n K II n V II II 11 1< IS M II n i> f •' 'i <i » is u <i n <i si si si si st ss si si st s> ii ii ii ii it is ii i> u ii ii ii ii ii it is ii ii ii ii n 

SIECC-9(Mi 


NAME 


COLUMN FORMAT DESCRIPTION 


CULL 


Used to delete individual 
measurements from the solution 


11-25 


E15.6 


The number* of the first measure- 
ment in a series to be deleted. 


26-40 E15.6 The number* of the last measure- 

ment in a series to be deleted. 
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CULL (cont.) 


IF CARD OMITTED: No measurements are deleted. 

NOTE: A maximum of 50 series of points per arc may be 

culled. The series may fall in any order. Only one 
series per card is permitted. If only one measurement 
is to be culled, then that number must appear both as 
the start and stop number of the series. 


*The numbers which correspond to each measurement are printed 
to the far left on the printout. 





DATA 


1)ATA 


1 1 II 


II 

I ! ) 4 S I 

I 1 1 


Hi 


m2 2 2|2 

3 3 3 3 3 
4 4 4 4 4 
S 3 3 S 3 S 

ilMS 

n } ni\ 

1 1 1 1 1 


2.1 


ioiiiiitiiiiiiiiiiiiiiiiiiiiioiii<iiD>iiiii*imoi************ ********** 

I I II II \1 U 14 li II '* ! !l If TI 71 7J 3* B H 3J II II M ?’ 31 « I* li •* 31 3' *' H *3 M li 4E 4* (» IJ II SI I! SJ W il SI 51 il Si ■ II II II *4 B ■ II II ■ H II H IS 14 M II I 

n 1 1 n 1 1 n 1 1 1 n 1 1 1 1 1 1 1 1 I M 1 1 It n 1 n 1 1 1 1 1 1 1 1 1 1 1 1 II 1 1 11 1 1 11 1 1 M 1 1 n 11 


3 3 3 


2 2 2|2 

3 3 3 

4 4 4 4|4 

5 5 5 
III E E 

7 7 7 711 

III I I IM 


22222222222222222222222222222222222722222222272222222222222222 

3333333333333333333333333333333333333333333333333333333333333 


55555555535555555555555555535555555555555555555555555555555555555 

E I E E E E E 6 E 1 1 1 1 « fi fi E E E 5 1 1 I fi S E E I fi E fi fi fi S E fi E E fi E E E fi E S E E fi E ; fi e E E E E E E E 1 1 E S 6 E 
7 7 7 7 I 7 7 7 7 7 7 7 : 7 7 7 7 7 7 M 7 7 7 1 7 7 7 7 7 I 7 1 1 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 

I 1 1 I 1 1 1 1 1 1 I 1 1 1 1 1 E 1 1 I I I I II I I I I I I I I I I I I I I I I I I I I I I ! 1 1 I 1 1 1 1 1 1 1 I 1 1 1 


nillStllSli9l99lll!IISMIIlll9HI!l''9lll9llllllll*************************** 

* *1 l| ’ I » l|(ii II i] 14 IS U '1 II ’I II 71 27 7? li II .*1 71 If 3*“ ’’ .1 -J ■* 3^ I*! I' ’’ *' ** ** •' • ' '* 5? SI 5* 5S SI SI 'A S3 Sfi !• I.’ •«. |l IS H 1^ II 63 18 D >2 *3 i* IS ’I 

•necc-boa* 


NAME COLUMN FORMAT DESCRIPTION 


DATA 1-4 A4 Specifies type of input data 

tape . 

11-25 E15.6 Integer part of number indicates 

type of tape used = INDATA: 

1 = residual tape 

2 = ARC tape 

Decimal part of number indicates 
file of ARC tape on which arc 
is found = IFILE. 


IF CARD OMITTED: Default value of INDATA = 1 is used. 

IFILE is not used. 


nil 

I n n K 
I I I 7 

! 2 2 2 

13 3 3 

14 4 4 

15 5 5 

IEEE 

'777 

I E I I 
19 9 9 

• I* >1 ic 
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EDIT 



cflIT 


ODD soDOOoooodiiooDoiBooonoooiieiensooDBeioiQaooooooagoooooiiooDQoaQDooooooogoooo 
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1 I ^ECC-50a< 


NAME COLUMN FORMAT DESCRIPTION 

EDIT 1-4 A4 Determines editing threshold 

level . 

11-25 E15.6 Edit threshold value in meters. 

IF CARD OMITTED: Default value of EDIT = 50.0 is used. 
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I 


PLOT 


000 OlOOe OOOOOOOOODIieOODOOOIOIBOIOOOOIOOOOOOOOOOOOOOOOOOOlOQODOflllOOOOOIIOeOlO 
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3 5 S S 5 S 5 S S S 5 S 5 5 S 5 5 S b 3 5 S S 5 S 3 S S S 3 5 5 S S 3 S S 5 5 5 S 5 5 5 S 5 5 5 5 5 5 5 5 5 S 5 S 5 3 S 5 S 3 3 S 5 S 5 3 5 9 S 3 5 S 5 S 5 S 3 

to 0 I fi 6 S S i S 6 E E 0 ( C 6 G 6 6 6 0 t S E 8 E S B S E E ( t C B 6 B S ( 0 0 S 1 1 0 6 S 6 G 6 6 0 E 6 6 6 6 E 0 S i 0 0 ( t i S S ( t 0 S 0 0 t fi t E 

11 in m m 11 1 1 ij ni n : 1 1 11 1 in 111 n n n 1 m m ni ni n 1 m 11 n 1 1 11 111 n 111 n 11 

1 1 1 1 1 I I S 1 1 I I 0 I E 0 8 0 I 0 I 8 0 t 0 I 8 B 0 1 1 1 I I 8 8 8 0 0 0 I 0 1 1 I B B B 8 B B 8 0 1 1 1 0 I t ! 8 1 1 I 8 1 1 1 1 1 1 1 8 I 1 1 1 1 0 E 

9 9 9 9 9 9 9 9 9 ! 9 8 9 9 9 8 9 9 9 n 9 S 9 9 8 9 3 I 9 9 9 9 S 9 9 ! 9 8 a 3 ! 8 8 8 9 S 9 9 9 S 9 S 9 9 0 9 9 9 9 9 9 9 9 9 8 9 8 S 8 9 9 9 S 9 9 8 9 ! 9 

I ) I I s 1 ' I s I mm ms im II ii ii ji n ii ii is ii it is is j’ c ii ii i> is ii ii ii i< i' < •: n <i <s ii i' <i is si s’ si si si ss m si s> si is ii ii ii ii ts it i' ii is io n n ii n is ii ii ii is ii 

MECC-SOSI 


NAME 


COLUMN 


FORMAT 


DESCRIPTION 


PLOT 


Determines if output tape is to 
be generated containing variables 
time, measurement, forward filter 
results, backward filter results, 
smoothed results, residuals of 
measurements about smoothed 
results. 


5 


PLOT (cont.) 


NAME COLUMN 


7-10 


IF CARD OMITTED: 


^Logical units 1, 
units . 


FORMAT DESCRIPTION 

14 I£ greater than 0, indicates 

unit* to be used for writing 
tape = IPLOT. 

Default value of IPLOT = 0 is used and 
no output tape is generated. 

2, 11, and 12 cannot be used as output 


1 


I 


REVNUM 


«<EVNlM 1710. 73. 

II I llisiaooo 0 0 1 1 0 0 0 0 0 0 I 0 0 I I 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 M I B ■ M I • i l> I 0 0 0 fl I 8 

1 I II 5 I i I I u II inj « IS II II II II n !i n I! mi » II 2 » a 31 II a « H IS a ); «a « <Mi < 1 M « mi n 4IH SI « a n 5s * SI SI SI mt B u II B ■ II ■ « II II ri )] H is ii ii ii ii ■ 

) n 11 1 1 1 1 1 1 1 I n 1 n 1 11 11 1 1 H 1 1 1 1 1 1 1 1 1 1 n I n 11 n M n 1 1 11 1 1 1 11 1 1 1 1 1 11 11 1 1 1 1 1 1 1 1 1 1 
l 222222n222222222222222222n22?2J22222222212222m222222!2J22222222222n2222222 
213333333333333 33333333333 333333333333333333333333333333333333333333333333333 

4444 44444444444444444444444444444444444444444444444444444444444444444444444444 

5 S S S 3 S S S 5 3 3 5 5 3 S b S S 3 5 3 S S S S S 3 5 S S S S 5 5 5 S 3 S 5 S S S S 5 5 5 S 5 5 5 5 5 5 3 5 5 5 3 5 3 S 5 5 S S 5 S 3 5 S 5 5 5 5 S S S 
8 S 8 I t i t B i i I 8 8 B i 6 5 B S t 8 8 G S I ( 8 i 3 « 8 6 S E t 6 I 8 S S i i ( E E i ( B t B 6 6 6 6 t 8 S E 6 I 6 S i 8 E ( I 8 ( i ( 8 8 S E t E I ( ( 

] m m m n 1 n n n n n : 1 1 n i m n 1 1 n n m ii m n n n j ni in n i ii m in i n ii ii , 

I 1BIII1I1M8I888 818888188888 181111881111111818181111111111811111111111111111118 j 
Ml IS MinSSS 81881 1 SS8IISS 11 llSSSSlIllllllinilMinSSl! 8 9 8 1 3 8 9 9 8 8 8 8 8 8 8 1 1 1 1 1 9 8 1 

t 7 3 4 5 i ' I 1 M •! 12 I) ]« I) |( |] im n 71 JJ 73 U 23 » 7MI 71 V V 33 X 33 17 )f 19 if It «7 41 44 4) « 41 4| 4« M SI 57 S3 S4 S3 » S3 S( H H (1 (7 (3 «4 IS U S3 H II M 17 13 ]4 3S H 13 M 7S N 

^CC-SMi 


NAME COLUMN FORMAT DESCRIPTION 

REVNUM 1-6 A6 Inputs rev number and segment 

number to be used on SEAHT tape. 

11-25 E15.6 Rev number = REVNUM 

26-40 E15.6 Segment number = SEGNUM 

IF CARD OMITTED: Default values of REVNUM = 0 and 

SEGNUM = 0 are used. 

NOTE: If SEAHT card is present and ISEAHT is greater than 

0 then this card must be present to give correct 
values to REVNUM and SEGNUM. 
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I 


SEAHT 


0 DO 0 Q 0 0 0 Q 0 0 0 D 0 a 0 0 O 0 S 0 D 0 0 0 0 0 0 Q 0 0 1) D t 0 0 D 0 D 0 0 0 0 0 0 D 0 D D 0 0 0 0 0 0 0 0 g 0 0 0 0 0 0 C OJ 0 0 0 II D 0 0 0 0 0 U 

I i : I i t ' I 9 II II I! i3 14 IS l( II It I) 21 21 22 2] 21 29 21 21 21 29 34 31 32 33 34 IS 3S 31 34 39 II 41 42 43 41 IS 4( 41 (I 49 S9 SI S3 33 SI SS SE S2 SS SS 13 (I 12 (3 El CS 51 51 51 51 39 II 32 33 II .’3 |E P It II 43 

II I n 1 M 1 11 1 n 1 11 n n M 1 11 11 11 n I 1 1 n n 11 I n 1 1 1 1 1 1 1 11 1 n 1 I t i I n n i I1 1 1 11 1 1 1 1 1 1 

nil 111 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ? 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 Z 2 2 2 2 2 2 2 2 2 2 2 2 2 2 

3333 333333333333333333333333333333333333333333333333333333333333333333333333333 


5 555555555555555535555555555555555555555555555555555555555555555555555555555555 

6 6 G 6 G 5 6 6 6 6 6 E 6 B 6 E E 6 E E 6 6 6 6 B G 6 B 5 D E 6 E 6 6 B 6 G 6 5 G E 6 6 6 E 6 6 G G 6 6 6 6 6 6 6 E g B 5 G 6 6 6 E E 6 e E 6 E 6 E E 6 G 6 B 6 

1 n 11 n m m 1 1 11 111 111 : n m 1 m ] m n 11 1 m m 11 11 n n m 1 n n 1 n n n n m n n 

ESS B 8 8 B B 1 8 G G B 8 E S a B e 8 8 8 B 8 8 B e 8 8 a E 8 8 8 B 8 8 B B e e S G 6 8 6 8 e 8 e 8 8 8 E S S e B S B 8 B B 8 8 B e 8 B S 8 e 8 8 8 8 B B 8 

99939998999989933903999899999999993999999999993999999999339399999999993989399399 
I 4 < j > » I •’ '1 13 II ir II I 'I 1 :a 31 j: :■> ’1 fi 31 :i ir n 33 31 11 is is 11 11 1' if ' 'i 11 « is n i' ‘i n n s si si si si ss s' ;s 59 cj 11 ii 11 ti is ti s' 11 ss n 11 'i 11 1« is >i n u '9 ii 

NECC-S08I 


NAME 


COLUMN 


FORMAT 


DESCRIPTION 


SEAHT 


Determines i£ output tape is to 
be generated for use as input 
to SEAHT. 


7-10 


If greater than 0; indicates 
unit* to be used for writing 
tape = I SEAHT. 


IF CARD OMITTED: Default value of ISEAHT = 0 is used and 

no output tape is generated. 


Logical units 1, 2, 11, and 12 cannot be used as output 
units. 
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SIGMA 


SlWfk 


9 1 0 a g| 

I i S < S I 
1111 1 


9 a 9 9 9|9 

I ) s < s 


II 0 0 9 0 0 0 9 C 0 0 0 0 0 I 0 lU 

' I t II II 12 13 14 IS II I) II II 21 21 22 23 24 2s|l 

1 I I 11 I 1 1 1 1 1 11 I 1 1 1 1 1 


9 9 9|9 

I ! 


2.00 


2 2 2 2 2|2 2 2 2|Z 2 2 2 2 2 2 2 2 2 2 2 2 2 2 
333333333333 3333333333333 
444 4444444444444444444444 
S a S S S 5 5 5 5 5 5 S 5 5 5 S S 5 5 5 5 S 5 S S S 

{ « I s s a 1 6 i 9 9 1 1 a ( 1 1 1 i ( I s a 1 6|t 
n m m m ni n n n n :i 

1 1 1 1 1 I I M 1 1 8 19 18 1118 19 8 


1.00 


0 0 0 1 1 1 1 0 1 1 0 I 8 I 0 1 1 1 1 1 • M I 0 0 0 0 I 0 8 0 0 0 0 0 0 I 0 0 0 0 0 0 0 8 0 8 0 0 I 0 

21 21 21 2i 31 31 32 a II a ■ II M M 41 41 42 41 44 41 4i 41 41 41 It SI 12 S3 S4 SS Si SI SI SI H II 12 13 14 IS 14 II II II II II 12 |2 14 IS II II N H a 

1 I 1 I 11 1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 1 1 I 1 11 I I 1 1 1 1 1 11 1 1 1 1 1 1 1 1 1 1 1 1 


22222222222 
3 333333333 

44444444444 

5 s s 5 s s a s a s s 

9 9 S « 9 9 9 8 9 9 I 

In 1 7 7 1 7 7 1 7 7 7 

888088888 




99999999999 


2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 Z 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 
333333333333333333333333333333333333333333 

444444444444444444444444444444444444444444 

a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a a 

966696899991699999999999999999999999999999 

777177777777777777777777777777777777777777 

8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 ! 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 


999999999999999399999999999999999999899999 


9 9 9 9 9 9 9 9 9 9 9 9 9 9|9 

l> I! |] 14 IS II II II II II II II II II IS|II II II II Il> II If II II IS » II « 14 4. <I 41 44 IS II II 41 H H SI St SI S4 SS II SI Si SI H II II <1 II IS tl II II II II II II II II IS II II 14 II II 
4IECC-5C»f 


NAME 


COLUMN 


FOW4AT 


DESCRIPTION 


SIGMA 


1-5 


A5 


Determines (i) Geoid undulation 
sigma and (ii) Altimeter noise 
sigma. 


11-25 E15.6 Geoid undulation sigma in meters 

SIGEO. 

26-40 E15.6 Altimeter noise sigma in 

meters = SIGNO. 


IF CARD OMITTED 


Default values of SIGEO = 2.00 and 
SIGNO = 0.60 are Used, 
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1 i 


8.3 SPECIFIC EXAMPLES 

8.3.1 Sample Deck Set-up and Optional Input 

The following is a sample of a working deck set-up. 
This run would operate on sea surface heights obtained from 
file 1 of tape 5344 which contains the pass for REV 1710, 
segment number 73. This job will produce two output tapes. 
The output tape used for plotting is written on unit 20 and 
the tape used as input to SEAHT is written on unit 21. The 
default values for CORLEN = correlation length, EDIT, 

SIGEO = geoid undulation sigma, and SIGNO = altimeter noise 
sigma are used because of the absence of their respective 
keyword option cards. Also, no measurements are culled. 


I 


COL. 


1 


•>> 

i 


I 


$ 

$ 


8 11 
IDENT 
OPTION 


1 ^ 26 
300902, REV1710 
FORTRAN 


PROGRAM BINARY DECK 


$ 

EXECUTE 


$ 

•LIMITS 

20,30K, ,10K 

$ 

TAPE 

02,X2D, ,5344 

$ 

TAPE 

20,X3D, , , ,SAVE 

$ 

TAPE 

21,X4D, , , ,SSHT 

$ 

FFILE 

02,NSTDLB,MLTFIL 

$ 

FILE 

11,F6S,20R 

$ 

FILE 

12,F7S,20R 

$ 

INCODE 

IBMF 

PLOT 

20 


SEAHT 

21 


REVNUM 

1710. 

73. 

DATA 

2.1 


$ 

ENDJOB 
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I 


8.3.2 Sample Output 


t 

I 


Specific examples of the program output are reproduced 
on the following pages. A sample of the optional input parameters 
and a portion of the program output are reproduced in Tables 
8.1 and 8.2. The RMS can be found printed out after the processed 
results. Table 8.3 is a list of the indices corresponding to 
the places where the state transition matrix is recomputed. 

These indices indicate where time gapes occur in the data. 

Notice that the indices of the state transition matrix occur 
in pairs. The actual time gap in the data occurs between the 
first of these data points and the previous data point (e.g., 
in table 8.3 a time gap occurs between data points 1717 and 1718} . 
Table 8.4 is an example of the output produced when data points 
are culled. In addition to the list, points which are culled 
are tagged with an asterisk on the printout for each recognition, 
as shown in Table 8.5. Equivalent messages and tags are also 
output when data points are edited (i.e., deleted). 
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edit* 5Q,00 

GEQID UNPUUATIOr:) §IGMAb ?,o? 
AUTIHETER NOISE S|,GMA* ?|50 
CORKEuATION WfcNGTMs 5>g,09 


as 


Table 8.1 


(NHX 

Tl*^E OF 

DaTa 

Elapsed 

TjHE 

KEASU8MENT 

FOftWARO 
HE IGMT 


xy^MoD 


ss, 

ssss 

CSEC> 

(METERS) 

t meters ) 

?1 

75o7o6 

1638 

36 

1 6657 

65,8465 

-45,6250 

-45,3593 

?2 

75p7o6 

1638 

36 

,9684 

65,9489 

-46,2344 

-45,3121 

?3 

75p7p6 

1638 

J7 

, 070=^ 

66, 0513 

"46, 4453 

-45, 3001 

?4 

75p7|35 

1638 

37 

|173? 

66,1537 

-46,9219 

-45,3599 

?5 

75o7q6 

1638 

37 

,2756 

66.2561 

-47,6094 

-45,5l69 

V6 

750^05 

1638 

37 

,3780 

66,3585 

-47,4375 

-4?, 6368 

V 

75o7p6 

1636 

37 

,4804 

66,4609 

-46,9141 

-45 ,6700 

?8 

75c7p6 

1638 

37 

|5825 

66,5633 

-45,B5l6 

-45,5436 


75o7p6 

1638 

37 

,6852 

66,6657 

-46, 7578 

-45,5594 

190 

75p7o8 

1638 

37 

, 787 6 

66,768i 

-46,6464 

-«5,557B 

101 

75p7p6 

1638 

37 

|0’o: 

66.87q5 

-46,6797 

-4S,56j& 

102 

75p7p6 

1638 

37 

.9924 

66.9729 

-47, UOOO 

-45,6101 

103 

750/06 

1638 

36 

1 0948 

67,0753 

-47,4922 

-45,7276 

104 

75o/n6 

1638 

30 

|197h 

67,1777 

67,2802 

-47,4922 

-45,8356 

185 

75o7p6 

1638 

38 

,2996 

•46 , 46q9 

-45,7826 

106 

75q7o6 

1638 

38 

|4c20 

67,3826 

67.4850 

-45,9375 

-45,6550 

107 

75(|/(36 

1638 

38 

,5o44 

-47 , 14o6 

-4»,7132 

108 

75o / ()6 

1638 

38 

,6068 

67,5874 

-4b, 7 q31 

-45,7021 

109, 

750706 

1638 

38 

,7092 

67,6898 

-46,1094 

-45,6036 

110 

75p7p6 

1638 

30 

r6n7 

67,7922 

"46,0156 

-45,4989 

111 

750706 

1638 

38 

|914i 

67.8946 

-46, 789i 

-45,5160 

112 

7507Q6 

1638 

39 

,0165 

67,9970 

-46,1406 

-<5,437o 

113 

75fl7o6 

1638 

39 

}1187 

66,0994 

-46,6594 

-4>,4337 

-45,4282 

il< 

75p7(36 

1636 

39 

»i2l8 

68,2018 

-46,5859 

115 

7 5 0 7 C 6 

1636 

39 

,3237 

68,3042 

-47,2613 

-45.5264 


75o7p6 

163» 

39 

,4261 

68 ,4066 

-46,0625 

-45,4386 

117 

75o7o6 

1638 

39 

,5265 

6B',5o9o 

-46,5938 

-45,4363 

118 

75o7q6 

1638 

39 

,63q7 

68,6114 

-46,8047 

-45,4661 

1*9 

75q7o6 

1638 

39 

,7333 

68,7138 

68,8162 

-46, 89o6 

-45,5d7o 

-4J,5lf3 

120 

75p7o6 

1638 

39 

|0357 

-47,2578 


BACKWARD 

smoothed 

residual 

MEIQHT 

HEIGHT 


jmetersj 

cmetehs) 

(METERS) 

-45,5060 

-46,6794 

1.0544 

-45,67i5 

-46,6645 

0i4501 

•45,7529 

-465.6890 

0.2437 

“45 , 8 064 

-46,6930 

•0,2289 

-45,7873 

-46,6967 

•0,9127 

-45,6564 

-46,7001 

H0.7374 

-45,5420 

-46,7034 

•0,2107 

•45,5o20 

-46,7065 

0.8549 

-45,6282 

-46, 7q95 

*0.0483 

•45,6205 

-46,7125 

0,0640 

-45,6297 

“46,7155 

0,0358 

-45,6346 

-46, 7i86 

•0.26l4 

-45,5896 

-46,7219 

r0.7703 

-45,4633 

-46,7255 

•0.7667 

-45,3263 

-46,7294 

0.2684 

-45, 3487 

-46, 7336 

0.7961 

-45,4560 

-46,7382 

•0. <025 

-45,3828 

-46,7432 

0.0400 

•45,3754 

-46,7485 

0.6392 

-45,4640 

-46,7543 

0.7387 

-45,5765 

-46, 7605 

r0.0286 

-45 ,5764 

-46 , 767o 

0.6264 

• 45 , 68q 4 

•46,7736 

0,1644 

-45,7i69 

-46 , 78o9 

0.1’50 

-45,7645 

-46, 7883 

•0. 4930 

-45,7033 

-46,7958 

0,7333 

-45,8310 

-46,8033 

0.2095 

-45,8840 

-46 , Bio7 

0.0060 

•45,9o6B 

-46,8160 

•0.Q727 

•45, *166 

•46,6249 

•0.4329 


SLOPE 

SIQHA 

SIGMA 


SLOPE 

SMOOTHED 

(ARCSEC) 

(ARCSEC) 

(ME7EHS) 

-1.3416 

3.3574 

0.1162 

•1.2151 

3,3574 

0.1182 

•1.1055 

3.3574 

0.1182 

•t.0086 

3.3574 

0.1162 

-0.9222 

3.3574 

0.1182 

-0,6469 

3.3574 

0.1182 

•0,7841 

3.3574 

■‘0.1182 

•0.7358 

3.3574 

0.1182 

-0.7033 

3,3374 

0.1182 

-0.6877 

3.3574 

0.1182 

-Cl 69 q3 

3,3574 

0.1182 

-0.7132 

3.3574 

0.1182 

-0,7584 

3,3574 

0.1182 

-0.8283 

3.3574 

0.1182 

-0,9246 

3,3574 

0.1182 

•1,0471 

3,3574 

0.1182 

•ltl936 

3.’ 3574 

0.1182 

-1,3599 

3,3574 

0.1182 

•1.5408 

3,3574 

0,1182 

•1,7293 

3,3574 

0,1182 

-1.9172 

3.3574 

0.1182 

-2,0960 

3.3574 

0.1182 

-2,2575 

3.3574 

0.1182 

-2,3938 

3,3574 

0.1182 

-2,4984 

3.3574 

0.1182 

•2,5655 

3.3574 

0.1182 

-2.59Q9 

_ 3:357; 

0.1182 

•2,57ie 

3.357< 

0.1182 

r2,5042 

3.3574 

0.1|82 

•2,39o7 

3.3574 

0.1182 


e 8.2 


SyUMB = 49136, ACTIVIT¥ # • P2o RtPQRT CODE 

state transition matrix f 1 
state transition matrix ■ 2 

bUTfc transition MAtRiX f 364 
STaTE TRA^J5IT10^' MATRIX ? 365 

glATt transition matrix f 781 
state transition matrix » 732 

state transition matrix f 1322 
SUTe transition matrix f 1323 
STaTe transition matrix f 1718 
state Transition matrix » 1719 
SUTE transition matrix ? 1994 

state transition matrix » 1999 
SUTE transition matrix f 2264 
state transition MAT''*1X f 2265 
STATt Tra»;SITION matrix * 2326 

SUTe transition matrix • 2327 

state transition matrix 9 2409 

state transition matrix » 2410 

O, state transition matrix f 2436 
O' state transition matrix • 2437 
SUTe transition matrix » 298I 
state transition matrix * 2882 
STaTE transition matrix f 3098 
SUTE transition matrix ■ ?099 

STaTe transition matrix • 4243 

state transition matrix 9 4242 

STaTe transition MaTRIX i 3097 
SUTe transition MATRIX 9 3q96 

SUTe transition matrix » 2880 
state transition matrix f 2079 
SUTe transition matrix 9 2435 

STlft transition matrix f 2434 

SlATt transition matrix f 2408 

SUTe transition matrix « 2407 

SIate transition matrix « 2325 

STATt transition MATRIX 9 2324 

state transition matrix f 2263 
ST.AtE transition MATRIX 9 2262 

St.AfE transition matrix » 1993 

SIaTc transition matrix 9 1992 

state transition matrix f 1717 
state transition matrix f 1716 
state transition matrix ■ 1321 


^2i RECORD COUNT ■ 00004B 



Table 8.3 


SNUHa « 49592, activity $ * 02# RcP0«T CODE » 06? RECORD COUNT « 005459 


MEASUREfiENT NO, 261lCULl£D 
measurement no. 2612CULLED 
measurement no, 2613CULLED 
measurement no. 2614CULLED 
measurement no, 2615CULLED 
measurement no. 2840CULLED 
measurement no, 2B41CULLED 
measurement no. 2842CULLED 
measurement no, 2843CULLEP 
measurement no. 2644CULLED 
MEASUREMENT NO. 2B45CULLED 
measurement no, 3180CuLLED 
measurement NO. 3181CULLED 
measurement no, 3182CULLED 
measurement no, 3183CULLED 
measurement no, 3154CULLED 
measurement no, 3504CULLED 
measurement NO. 3505CULLED 
measurement no, 3506CULLED 
measurement NO. 3507CULLED 
measurement no. 3508CULI,ED 


INDX 

time 

or 

DATA 

elapsed 

time 


yymhdd hhmm 

SS.SSSS 

ISEC) 

2P21 

75o93l 

58 

24.9293 

422.9346 

2822 

75o83l 

58 

25.0318 

423, o37o 

2823 

75o83l 

58 

25.1342 

423.1394 

2824 

75(3031 

58 

25;2366 

423,2418 

2825 

75o93l 

59 

25.3390 

423.3442 

2826 

75q33i 

58 

25.4414 

423.4466 

2827 

750831 

58 

25.5438 

423.5490 

2828 

750831 

58 

25.6462 

423.6514 

2629 

75o83l 

59 

2,5.7466 

§23.7538 

2830 

75q83i 

58 

25;8510 

423.8563 

2831 

750831 

58 

25.9534 

§23,9587 

2832 

75o83l 

59 

26.0558 

§24 . o6ll 

2833 

75o83l 

59 

26,1582 

§24.1635 

2834 

75o33l 

58 

26.2606 

424.2659 

2835 

750031 

59 

26.3630 

424.3603 

2836 

75q831 

58 

26.4654 

424 . A7q7 

2837 

750831 

59 

26:5676 

§24,5731 

2838 

750631 

58 

26.6702 

§24.6755 

2839 

750631 

59 

26.7726 

424.7779 

2840 

750831 

68 

26:6751 

424.8603 

2641 

750831 

59 

26:9775 

*24.9827 

2842 

750831 

58 

27.0799 

425,0851 

2643 

75o83l 

59 

27; 18 23 

425.1875 

2844 

750831 

59 

27.2847 

§25;2899 

2885 

750831 

59 

27,3871 

425.3923 

2846 

75o83i 

58 

27.4895 

«25,4947 

284 7 

750831 

59 

27:5919 

§25,5971 

2848 

750831 

59 

27:6943 

§25.6996 

2849 

750831 

58 

27,7967 

425.6020 

2850 

75o83i 

59 

27:8991 

925.9044 


MEaSURMENT 

forward 

height 

backward 

HEIGHT 

Smoothed 

HEIGHT 

fMETERSJ 

<METERS) 

tMETERS) 

(METERS) 

-i5.25oo 

= l4'.4io3 

-16,3058 

•19,6341 

-15,1563 

514.4837 

-16,5i7i 

tl5.7oi7 

-14.3594 

814.4339 

-16,7538 

’15,7646 

-l5.o39i 

514,4878 

-1T;1276 

’15.8221 

-19,8672 

=14,6594 

-l7,*ll9 

-15,8738 

-16.7578 

•14.9380 

-17,5745 

-15.9i93 

-19,4375 

=i5,59iq 

-17.5943 

’15.9584 

-20.57o3 

»i6,354g 

-I7,l7i5 

-15.99o9 

- 21.8281 

=17,2339 

-i6;5i35 

-15,5728 

’16.0171 

-20.9219 

=17.8996 

ni 6 .0369 

-l8,l0l6 

il8.o82i 

-l4,674i 

’l6 , o5o5 

-14.3359 

=17,6777 

-l4,i484 

’l6,0580 

-13,1641 

8lT.ll4g 

-14.2235 

-16.0998 

•13,4922 

=i8.627i 

-l4 ,53 o8 

-16.0962 

•14,8438 

=16.3656 

-14,82Bi 

■:16.0478 

-16.8399 

“l8.4092 

-14.8947 

-l 6 , 0353 

-15.9688 

=16,3i98 

-14,4972 

-I®. 0194 

•i4,1641 

=15.9589 

-14,i946 

- 16,0008 

•14.3047 

=l5.6466 

-14,3233 

’15,9799 

-14.9638* 

-15,5386 

-14,4452 

-15,9971 

-15.4141. 

P15.4224 

-14,555o 

•15,9330 

-l5 , 0078. 

•15.2987 

-14,6508 

-19,9080 

-l4 , 7Bi3. 

=l5.l68o 

-14;7522 

-14,8593 

-15,8629 

-15,3828. 

-i5.o307 

•15,8569 

-i 6 ) 3628* 

•14,0873 

-14,9519 

-15.83i7 

-l5.78i3 

•l5 , 0479 

-15,0399 

’l5,8o79 

-1§ . 67l9 

=14,8531 

-15,0058 

115,7846 

-14,9688 

•14.7686 

-15,1478 

•15,7639 

-15,4922 

•14.8018 

-15,2550 

•15,7449 

-15,2969 

• 14.-7846 

- 15,2908 

415.7287 


RESidUaL 

SLOPE 

SlRMA 

sigma 



. SLOPE 

smoothed 

meters') 

(ARCSEC) 

(ARCSeC) 

(METERS) 

0.3841 

-22.69§6 

3.4301 

0.1186 

0.5454 

t21.37i6 

3,4446 

0! 1169 

1.4052 

-19,7644 

3.4599 

0.1192 

0.7830 

-l7,945i 

3.4756 

0.1196 

0 .0066 

•15.8756 

3.4911 

0.1201 

’0.8385 

-13.60*6 

3.5058 

0.1207 

r3.479i 

-11.1933. 

3^5190 

0.1214 

’4.5794 

-8,6859 

3.5300 

0.1222 

45,0110 

-6, 1568 

3.5380 

0.1231 

n§ ,885o 

-3.6803 

3.5424 

0»12*1 

42.o5ii 

-1.3272 

3.5427 

0‘1252 

1.7221 

0,8429 

3.5383 

0*1263 

2;B957 

2'.7B74 

3.5292 

0.1275 

2,5640 

4.48i5 

3,5i54 

0*1287 

l,204o 

5.9i35 

3.4975 

0.1298 

r0.6006 

7, 0855 

3.4763 

0-1309 

0.o5o7 

7,9999 

3.4530 

0-1319 

1.8367 

8.6692 

3.4294 

0.1327 

1'.6752 

*.i097 

3.4072 

0.1335 

0.9884 

9.3405 

3.3884 

0.134Q 

0.5190 

9.3603 

3.3747 

0.1344 

0,9002 

9,2456 

3,3675 

0.1345 

1,1012 

8.9513 

3.3675 

0-1345 

0.4741 

8.5113 

3r3747 

0.1344 

•0^5511 

7.9395 

3.3864 

0.1340 

0.0262 

7.2512 

3,4073 

0.1335 

1,1127 

6, 4640 

3.4294 

0.1327 

0.7948 

5,5984. 

3,4531. 

0,1319 

0,2526 

4,6769 

3t4763 

0.1309 

0.4318 

3.7220 

3,4975 

0.1298 


Table 8.5 

Note the asterisks after measurements 2840 - 2845, 
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8.3.3 PLOT Tape Format 

This type o£ output tape contains information on the 
variables produced as output from the filter. The tapes 
are 800 BPI standard label tapes. The main use of these 


tapes is in plotting certain output variables versus time, 
as exemplified in Figures 8.1 and 8.2. The data record format 
for this tape is as follows: 

WORD 

FORMAT 

DESCRIPTION 

1 

R 

Time in seconds from the first data point 
of the arc. 

2 

R 

The input measurement as read from the 
input tape. 

3 

R 

The value* returned from filtering in the 
forward direction. 

4 

R 

The value* returned from filtering in the 
backward direction. 

5 

R 

The final combined smoothed value. 

6 

R 

The residual of the measurement about the 
smoothed value (W0RD(2) - WORD(5)). 


These values correspond to the a posteriori values returned 
from the filter after improvement based on the individual 
measurement. 
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8.3.4 SEAHT Tape Format 


This type of output tape is a tape for use as input 
to the SEAHT program. The tape is an 800 BPI, standard label 
tape. It contains one header record and then subsequent 
data records. The record formats are as follows: 


4 


SEAHT INPUT FORMAT 
HEADER RECORD 


WORD NO. 

TYPE 

DESCRIPTION 

1 

I 

Year and day of year 

2-3 

DP(FRAMTI) 

Frame time of day (seconds) 

4 

I 

Segment number - must be input 
to ALTKAL 

5 

I 

Rev number - must be input to 
ALTKAL 

6 

I (MTYPE) 

Measurement type (40 or 41) 

7-19 

I 

Zeroes (used for making all 


records same length) 
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SEAHT INPUT FORMAT 
DATA RECORD 


WORD NO. 

TYPE 

DESCRIPTION 

1 

R 

Time from first data point 
(minutes) 

2 

R(LAT) 

Latitude of first smoothed sea 
surface height in record 

3 

R(LON) 

Longitude of first smoothed sea 
surface height in record 

4-13 

R 

Smoothed sea surface heights 

14 

R(THITE) 

Tide height 

15 

R(TREF) 

Tropospheric refraction correction 


16 R(RAGCAV) Average AGC 


17-19 R Zeroes for future additions 
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APPENDIX A 

PROCESSING TWO SETS OF UNCORRELATED DIRECT MEASUREMENTS 


I 


A number o£ alternative algorithms for obtaining an 
optimum estimate from two sets of uiicorrelated estimates, 
or direct measurements, are presented in this Appendix. 


Given: Two sets 

of 

estimates 

and estimate error covariances 


s 

[Pal “ 

E([x^-x] [x^-x]^) 

[Xbl 

9 

[Pfa] = 

E([x^-x] [x^-x]^) 

E([x 

a-x][Xb-x]’’) = 

[0]^ — ^the estimates are uncor 


related. 

[x^] and [x^] may either be estimates, or direct 
measurements . 


Determine : 

New estimate and covariance 
Solution : 

It is well known and easily proven that the optimum 
estimate and its covariance are given by 

[i] ' [p] ([Pal'^LXg] + [Pbr^[X|,13 
-1 -1 

[p] = ([p^] + [Pi,l ) 
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I i 1 I f I J I ■ 

These expressions are not in the most convenient computational 
forms. We shall give several alternative forms and show the 
only matrix inversion required is ([p^] + 

Alternative Formulations ; 

[p] = ([Parh[PbJ * [Pali tPh)'^) 

= [PbH[Pbl*[Pal3'hPal = [PaK[Pa]*[Pb])’hpbl 
' ([Pbl^tPal' [Pa^^ f 

= [Pal- [Pal C[Pal*[Pbl) '^[Pal < — > Variance Equation 
= [Pa]-[KHPal 

[K] “ [PaJC[Pal + [Pbl) 4 — > Optimum Gain 
= [Pl [Pbl'^ 

[i] = [pmParhx^]^-[Pbrhxbn 

= [p] t([Pal‘^+[PbJ'^-IPbJ'^5 t*a)*tPbl'^t*b^' 

= [Xa]*lP]fPb'’hfXb]-[Xa]) 
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Applications : 


1. [x^] and [x^] may be considered as the forward 

and backward filtered estimates. 


2. Application to altimeter measurements 

*”x 

I 

In this case one may take [x, ] = 


and [p, ] = 


£00 
0 £ 0 

0 0 R 


£->•00 


Therefore , 


([PalMp,]) 


l/£ 0 

0 l/£ 

0 0 


(R+P. ) 

^33 


-1 


£-4-0 


K = [Pa]([Pa]^[Pbl) 


-1 




R'*‘P. 


33 


0 0 

0 0 

,0 0 


13 


^23 


33 
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[K]([x^]-[x^]) 




) 


[K][p^] = 




One sees these are exactly the same equations relating the 
a priori and a posteriori estimates given previously in 
Section 4.2. 
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APPENDIX B 

DERIVATION OF SMOOTHER DISCRETE TRANSFER FUNCTION 

In processing the mid-portion o£ a long data sequence,* 
the forward and backward filter gains approach their steady- 
state value.* The smoother becomes stationary and may oe 
represented by a transfer function. Since the transfer 
function gives additional insight to the smoother behavior, 
it is derived in this Appendix. 

From Equations (4.4) and (4.10) one obtains the follow- 
ing equation governing the "a posteriori" forward state 
estimate , 

= [4.^1 [X^] - [K] (B.l) 

where 

[■(■^1 = ([E] - [KlolO]) [ 4 .]. 

[E] = identity matrix 

In the steady-state [4^] and [K] are constants and one may 
take the z- transform of Equation (B.l) to obtain 

z[X^] = [X^l + Z[K] , 

or 

f f ” ^ 

[X*] = z(z[E] - [4^) [K] (B.2) 

For constant data sampling rates and no missing or edited data. 
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I 


where 



I 


[X„] = Z{[x„n 

Similarly, the "a priori** backward state estimate 
is given by. 


- (♦) tx^,2] 

= [*] U**) [x^^j] * [K] (B.3) 

For brevity in the above equations we have suppressed the 
superscript **+’* for **a posteriori** estimate, but have 
retained the superscript **-** for **a priori*' estimate. 

By taking the z-transform of Equation (B.3), one 

obtains 

= [♦] 
or 

Zftx^:^]} = C(E] - '' [K] tB.4) 
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/ 

/ 

The smoothed estimate is taken as the optimally- 
weighted combination o£ the ”a posteriori" forward estimate 
and the "a priori" backward estimate; i.e., 

= [W]p ^ [W]g [xb;jl 

Taking the z- transform of this equation, and making use of 
Equation (B.l) and (B.4), one obtains 

z(XJ = Yjj z^{awIp(z[E]-[4.*D ^+[W]g['t'] ([El-zt*^]) ^))[K], 
or, the transfer matrix, 

(H(z)] 4 1 [X ] = (z[W]_tF(z)] * [W]g[F(i)])[K] CB.5) 

n 

where 

f 

F(z) = (z[E] - [4)^]) . 

The third component of [H(z)] is of particular interest as 
it is the z transfer function relating the geoid undulation 
estimate to the altimeter measurements. Notice that the 
transfer matrix contains negative as well as positive powers 
of z, indicating the state estimate depends on "past" as 
well as "future" data, which is, of course, the essence of 
smoothing. The gain [K] and weighting matrices [W]p and [W]^ 
appearing in the expression (B.5) for the transfer matrix 
are obtainable from the steady- state solution of the filter 
variance equation. 
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NOMENCLATURE 


= matrix 


= expectation of ( ) 


= altimeter noise correlation function 
decay constant, 1/b is the (1/e) noise 
correlation time 


= geoid undulation 

== converted measurement noise variance 
= optimum gain matrix 
= altimeter noise 

= estimate error covariance matrix 

= increment in covariance matrix resulting 
from dynamic noise excitation 

= altimeter noise variance 

= time 

= sampling period 

= white noise generating exponentially 
correlated altimeter noise 


geoid undulation state variables as de- 
fined in Equation (3,5) . x^ is the geoid 
undulation. 


T 


I 


1 


w 


altimeter measurement with orbit height 
and other known corrections subtracted 

white noise exciting the geoid undulation 
state 


a 


= 1/23 

= geoid undulation correlation decay constant, 
2.903/3 is the (1/e) geoid undulation cor- 
relation time 


|.V 


white measurement noise in converted 
measurement defined in Equation (6.4) 


[ 4 >] 


geoid undulation state transition matrix 
defined in Equation (3.7) 


= a priori geoid undulation variance 


Subscripts 


= time 

i,j,k = indices for time 

b = output from backward filtering 

s = output from smoother 


Superscripts 

+ 


immediately after measurement 
immediately before measurement 
estimated quantity 
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